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Summary 


The method of studying water movements by measuring the 
electrical potential produced by the transport of water through the 
Earth’s magnetic field has been applied to flow through the Irish Sea. 
Continuous readings of these potentials were made between 1954 
November and 1958 March on a submarine cable which lies on the sea 
bed of the North Channel between Portpatrick and Donaghadee. 
Similar measurements were made further south on a Nevin—Howth 
cable during the period 1956 December to 1958 March. ‘To calibrate 
the cable signals, direct current measurements were made at different 
depths in a number of localities in the North Channel. No measure- 
ments were carried out to calibrate the Nevin-Howth cable, however, 
but an estimate was made from values given in Admiralty publications. 

Daily values of the current in the North Channel referred to the 
long period mean as zero are presented and it is shown that high 
correlations existed during the winter months between these flows and 
the regional winds. During the summer, however, the correlation was 
not as strong. Linear regression equations relating the mean flow to the 
regional wind were computed for bi-monthly, six-monthly and total 
periods. The results showed that a wind from S 49° E was the most 
effective in producing a transport of water northwards through the 
North Channel and that under these conditions the flow amounted 
to about 1 per cent of the wind speed. Furthermore, estimates of the 
accumulated transport of water due to the wind, suggest that the mean 
flow over a long period deduced from the salinity distribution may be 
due almost entirely to the residual effect of the wind. From considera- 
tion of a current surge on 1955 November 2 it is shown that the flow 

lagged behind the generating wind by only about two hours. 
265 


e 
4 


F. Eviton, Inamam, RH. Ox & 
Earth’scrastin Antarctica and the surround- 


ton of polar geomagnetic 
The re- tization of a surface outcrop by currents. 
Small ia the theory of suriace — 


sper K. Samara & M. The potassium-argon age of the lava 
The eof the Skiddaw and Eskslale granites. 


of the orbit of an earth satellite. 
we E A new field stability. 
LEGG, A tlew test 

‘The palacomagnetic hist 

A Powis, J. & F. Nucleosynthesis during the early 
B. A. Manvssi, Eigenvibrations of the Egrth observed at Trieste. 


The ature of the Mohorovitic 
J A. Jacome, Oscillations of the 


following papers have been accepted far publication : 
mevement in Japan end ite vicinity. 
An analysis of the position of the Earth's magnetic pole in the 
A deviating permian pole from rocks in Northern Italy. 
Interaction betwee tide sad surge inthe Thames. 
Crane, Je, Heat flow in the Austrian Alps. : 
eer Lennon, ‘The deviation of the vertical at Bidston in response to the attraction Bins 
A. Savina & The of seismic 
t 


The Geophysical Journal 
of the 


ROYAL ASTRONOMICAL SOCIETY 


Vol. 5 No.4 October 1961 


The Flow of Water through the Irish Sea and its 
Relation to Wind 


K. F. Bowden and P. Hughes 
(Received 1961 January 4) 


Summary 


The method of studying water movements by measuring the 
electrical potential produced by the transport of water through the 
Earth’s magnetic field has been applied to flow through the Irish Sea. 
Continuous readings of these potentials were made between 1954 
November and 1958 March on a submarine cable which lies on the sea 
bed of the North Channel between Portpatrick and Donaghadee. 
Similar measurements were made further south on a Nevin—Howth 
cable during the period 1956 December to 1958 March. To calibrate 
the cable signals, direct current measurements were made at different 
depths in a number of localities in the North Channel. No measure- 
ments were carried out to calibrate the Nevin—Howth cable, however, 
but an estimate was made from values given in Admiralty publications. 

Daily values of the current in the North Channel referred to the 
long period mean as zero are presented and it is shown that high 
correlations existed during the winter months between these flows and 
the regional winds. During the summer, however, the correlation was 
not as strong. Linear regression equations relating the mean flow to the 
regional wind were computed for bi-monthly, six-monthly and total 
periods. The results showed that a wind from S49°E was the most 
effective in producing a transport of water northwards through the 
North Channel and that under these conditions the flow amounted 
to about 1 per cent of the wind speed. Furthermore, estimates of the 
accumulated transport of water due to the wind, suggest that the mean 
flow over a long period deduced from the salinity distribution may be 
due almost entirely to the residual effect of the wind. From considera- 
tion of a current surge on 1955 November 2 it is shown that the flow 
lagged behind the generating wind by only about two hours. 
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Correlation of the Portpatrick and Nevin records confirmed the 
occurrence of periods of wind induced flow through the Irish Sea as a 
whole. Uncertain changes in the zero of each cable, which were 


probably attributable to Earth currents, prevented a study being made 
of the water balance. 


1. Introduction 


It is known from the salinity distribution that, over a long period, there is a 
net flow of water through the Irish Sea from south to north. The mean velocity 
through the North Channel was estimated by Knudsen (1907), from considera- 
tions of continuity, to be 1-gkm/day. A more detailed calculation by the same 
method (Bowden 1950) indicated that a large part of the exchange of salt content 
through the North Channel was probably due to longitudinal mixing and that 
the mean flow was only about 0-5 km/day. This flow, small as it is, is much larger 
than the displacement current needed to remove the fresh water added to the sea 
by rainfall and run-off, which would be of the order of 0-02 km/day. 

Current measurements at mid-depth were made at the Skulmartin light-vessel, 
in the North Channel about 3-7km from the Irish coast, in 1935-36 using a 
Carruthers drift meter (Proudman 1939). Over a year the residual current was 
o-6km/day towards the south-east, suggesting that the flow at this point was not 
representative of that across the whole cross-section of the Channel. Individual 
monthly mean flows were sometimes north-westerly and at other times south- 
easterly, with maxima of 4km/day in either direction. Further measurements at 
the Skulmartin light-vessel, made during a 15-month period in 1953-4 with a 
Carruthers vertical log current meter, gave similar results.* There was only a 
poor correlation between the monthly residuals of current and wind. 

Other evidence that a southerly flow through the North Channel occurs near 
the Irish Coast, at least on some occasions, was provided by drift bottle results 
reported by Tait (1952). Drift bottles released within the Irish Sea in a number of 
investigations (e.g. Daniel & Lewis 1929) have shown that the drift of the surface 
water is variable and is determined mainly by the wind. The results have also 
given some indication, however, of a general northerly flow superimposed on 
the wind drift. Williamson (1956) discussed the value of various species of 
plankton as indicators of water movements in the Irish Sea and concluded that, 
in general, they confirmed the pattern of flow outlined above. 

The cable method of measuring flow through a channel makes use of the electro- 
motive force induced in the water by its movement through the Earth’s magnetic 
field. The potential difference across the channel is proportional to the mean 
velocity of the water across the whole cross-section. In 1953-54 a study was made 
of the flow of water through the Straits of Dover by this method (Bowden 1956). 
The recorded values of flow were correlated with the wind, which influenced the 
flow directly by its tangential stress on the surface of the water and indirectly by 
causing a piling up of water due to the constriction of the English Channel in the 
region of the Straits of Dover. The resultant slope of the sea surface gave rise to a 
component of flow. 

In the present investigation, the same method was applied to the flow through 
the North Channel, using a suitably located telephone cable. As in the case of the 
Straits of Dover, the voltage recorded on the cable was calibrated in terms of 
* Unpublished data, supplied by Cdr. A. L. Lawford, National Institute of Oceanography. 
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measured values of the tidal currents. To study the non-tidal flow, the tidal 
currents were eliminated and the residual flow across the section correlated with 
the wind over the area. The information obtained in such an investigation is of 
interest in connection with the exchange of water between the Irish Sea and the 
region outside the North Channel. 


2. Apparatus 


The observations in the North Channel were made on a telephone cable running 
from Portpatrick, Scotland, to Donaghadee, Northern Ireland, which was put 
at our disposal by the Post Office. The recording gear was housed in the telephone 
repeater station at Portpatrick. Records were obtained covering the period from 
1954 November to 1958 March. For a part of this period, observations were also 
made on a cable running from Nevin, North Wales, to Howth, Eire. The posi- 
tions of the cables are shown in Figure 1, and profiles of the sea bed across the 
channel at each location are shown in Figures 2(a) and 2(b). 


Fic. 1.—Irish Sea and the cable positions. 


The recording circuits at Portpatrick and Nevin. At Portpatrick the measure- 
ments were made on a reserve four-core submarine cable running from a point 
about two miles north of Portpatrick, Scotland, to Donaghadee, Northern Ireland. 
The circuit arrangements are shown in Figure 3. 

A phantom circuit was used as the return path for the water induced signal, 
the connections being taken from the centre point of the cable side of the line 
transformers associated with a side circuit. The sea earth referred to in the diagram 
is the unprotected cable sheath which acts as a continuous electrode across the 
channel. 
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A siphon-pen recording milliammeter, having a full scale deflection of 1mA 
and a nominal internal resistance of 1000 2 was used to measure the varying cable 
potentials. The instrument was of the centre zero type so that the ebb and flow 
of the tide were recorded as positive and negative movements of the pen across 
the paper. The chart speed was constant at 1 inch per hour, and this was adequate 
for the investigation. When the system was being set up resistances in the control 


(a) 


Distance (km) Howth - Nevin 
40 


(b) 
Fic. 2.—Profiles of the sea bed across the channel at each cable location. 


(a) North Channel (Portpatrick-Donaghadee cable). (b) St George’s 
Channel (Nevin—Howth cable). 


circuit were adjusted so that deflections at tidal springs would traverse most of 
the width of the recording paper. Twice a week a standardizing signal was applied 
to indicate an accurate zero and as a functional test of the recorder. The timing 
of the chart drive was also watched carefully and time marks were placed on the 
record to allow a comparison of the Portpatrick and Nevin records to be made. 

On the cable further south between Nevin in North Wales, and Howth in the 
Republic of Ireland, a similar set-up with slight modifications was used. The 
cables available did not terminate at the beach, but were brought straight to the 
repeater station at 1-2 miles from the shore, and as a result the sheath of the cable 
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was exposed to land earth connection for that distance. A block schematic diagram 
of the cable layout is shown in Figure 4. 

On the office side of the low pass filter at Nevin a blocking condenser, C, 
(8 uF) is fitted to prevent the d.c. flow from saturating the transformer windings. 
It was across this condenser, at A—B, that the recording device, similar to the 
one used at Portpatrick, was connected. As mentioned, a reserve submarine cable 
was used at Portpatrick, but at Nevin the cable was carrying telephone traffic 
during the time that recordings were being made. 


Portpatrick Donaghadee 
telephone repeater station telephone repeater station 


ale 


Local cable 


Fic. 3.—Recording circuit, Portpatrick—-Donaghadee cable. 


Howth telephone Submarine Nevin telephone 
cable 


repeater station 


Fic. 4.—Recording circuit, Nevin—Howth cable. 


3- Earth currents associated with magnetic disturbances 

For most of the time fluctuations of a variable nature, both in period and 
amplitude, were superimposed on the water induced cable signal. These extra 
variations were due to the earth currents associated with disturbances of the 
Earth’s magnetic field, and their occurrence has been discussed by one of the 
authors in a previous paper (Bowden 1956). It is interesting to notice, however, 
the widespread nature of these disturbances, which occurred regularly in the same 
form on both the Nevin—Howth and the Portpatrick-Donaghadee cables. A typical 
example is given in Figure 5, where the two records are shown side by side. ‘There 
is an obvious correspondence between the sets of fluctuations which resemble one 
another in detail. 


There were no recognizable time lags between similar magnetic disturbances 
recorded at Portpatrick and Nevin. When they were suspected, investigations 
showed that one or other of the recorder clocks were in error. This aspect, however, 
was not examined extensively, but the recognizable features always provided a 
check on the record positions. 

The Nevin cable is longer than the Portpatrick cable while the mean depth 
over the section is less, and this results in relatively stronger disturbances in the 
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Nevin to Howth cable. The fluctuations are, however, mainly of fairly short 
period and can easily be removed. Where there was any doubt about this elimina- 
tion the record was not included in the results. There were, of course, occasional 
magnetic storms when analysis of the trace was impossible (see Figure 6), but 
for most of the time the disturbances were small and a good record is shown for 
comparison. No attempt has been made to analyse the magnetic fiuctuations, 
they were merely removed from the record. 


4- Calibration of north channel cable records 


The calibration was carried out in terms of the tidal currents, by comparing 
the observed voltage on the cable record with measurements of the tidal currents 
made at the same time. Two sets of data were used: (a) Observations by H.M.S. 
Cook, at a position in the centre of the channel, during one complete tidal period 
on 1955 May 26; (b) Observations with drifting drogues, at three positions across 
the channel, using the M.V. William Herdman, in 1956 July. 


(a) Current measurements by H.M.S. Cook 


The observations were made with a 30 ft drifting pole logship, noting the speed 
and average direction of current over 30-minute intervals throughout the tidal 
cycle. The current reached a peak of 2-2kt (113.cm/s) in each direction. The 
cable record showed the voltage reaching its maximum values approximately 
30 min earlier than the observed current. There was no evidence at this time about 
the variation of current with depth or across the channel. Assuming the current 


values measured by H.M.S. Cook to represent the mean velocity across the whole 
section, the calibration factor would be 


1 volt = 3:23 kt (166cm/s). 


This gives an upper limit to the current, since the mean velocity across the whole 


section is almost certainly less than that in the upper 30ft in the centre of the 
channel. 


(b) Drifting drogues, observed from M.V. William Herdman 


Each drogue was in the form of two steel plates intersecting at right angles, 
of overall dimensions 3 ft x 2ft. The drogue was suspended by steel wire of dia- 
meter 0-075 in. from an elliptical buoy of horizontal diameter 14in and height 
1oin. Attached to the buoy was a bamboo cane 1oft high, carrying a flag of 
distinctive colour to mark it. Three drogues were normally released together, 
at depths of 5, 25 and sofm (9, 46 and 91-5 m) or 5, 50 and roofm (9, 91°5 and 
183 m). 

If the velocity of the water at the surface is different from that at the depth of 
the drogue, the drift of the buoy will differ from that of the drogue by an amount 
depending on the relative drag forces on the buoy, drogue and cable. From their 
dimensions, and assuming appropriate drag coefficients for the drogue, buoy and 
cable, the following correction formula was derived: 


Vp = tp—a(Vs—vp), 


where Vs = velocity measured by surface drogue (depth 5fm), vp = velocity 
measured by drogue at depth D, Vp = corrected velocity at depth D and « is a 


) 
| 4 
a 
7 
4 
vai 
j 


Fic. 5.—Corresponding disturbances on the two cables, 1957 April 26 
2000h to 1957 April 27 o8o00h. (a) Portpatrick, (b) Nevin. 


Fic. 6.—Comparison of two Portpatrick records to show: (a) a calm 
period with no interference due to fluctuations of the Earth’s magnetic 
field. (b) the effect of a magnetic storm. 
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(c) 1955 November 3. (Crown copyright; 
from the Daily Weather Reports issued by the Meteorological Office by 


permission of the Controller of H.M. Stationery Office). 
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coefficient, having the value 0-3 for a drogue at 25fm, 0°34 at 5ofm and 0-4 at 
100 fm. 

The buoys with drogues set at three depths were released close together, at a 
position fixed by taking bearings on landmarks with a sextant. The position of 
each buoy was determined subsequently, at intervals of 4 to 1 hour, by bringing 
the boat near to it and taking bearings in the same way. As long as the buoys 
remained in the general area of the station and within reasonable steaming distance 
of one another, they were allowed to drift freely. When they had drifted too far 
apart, usually after several hours, or had moved too far from the nominal position 
of the station, they were taken out of the water and released together again at the 
starting point. The main limitation of the method was that set by visibility and 
hours to daylight, since the William Herdman was not fitted with any radio 
navigational aid and all positions were fixed by visual observations. 

The positions of the stations and the periods of observation are shown in Table 
1. Owing to the limitations referred to above, it was not possible to obtain observa- 
tions over a complete tidal period in each position. However, at each of the three 
stations on the successive days 1956 July 9, 10 and 11, observations were made 
over the full duration of the ebb flow, and these three sets of data were combined 
to provide an estimate of the total flow across the section during the ebb. At the 
central station, No. 2, on July 25, the observations extended over g hours and 
included both the maximum flood and maximum ebb currents. 


Table 1 


Calibration of North Channel cable. Positions and periods of drogues observa- 
tions 


Station, depth, Mean position Date Period of Hours relative 
approx. position — — — July observations to HW 
Lat. N Long W hours 


Za 98m 54° §° 24°79 6 HW +4 to 
5 ml. NE of HW +6} 
Donaghadee 


2. 240m 54° 45'°6 HW —z to 
7ml. SW of HW +6 


Portpatrick 


130m 54° 47'°6 HW —z2 to 
4ml. SW of HW +6 
Portpatrick 


240m 54° 45°°6 HW —4} to 
7mi. SW of HW +44 


Portpatrick 


By plotting the successive positions of the drogues, the components of drift 
velocity were determined and corrected for the drag of the buoy and wire, as 
described above. The corrected velocity component along the channel, taken as 
positive in direction 150° true, will be denoted by U. The values of U at each 
depth, at a given station, were plotted as a function of time and the values at hourly 
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intervals relative to HW were read off. Since the three drogues were near the sur- 
face, approximately at mid-depth and near the bottom respectively, a simple 
arithmetic average was taken as the depth-mean U of the velocity component at 


each hour. The mean velocity perpendicular to a section across the whole channel 
was then taken as 


U = }(0, + 2U2+ Us) 

where 0;, U2, U3 denote the depth-means at stations 1, 2 and 3, the value at 
station 2 being given double weight because of the greater depth of water there. 
Observations taken on the consecutive days, July 9, 10 and 11, when the range of 
tide at Belfast was 9-5 ft, 9-4ft, and 9-3 ft respectively, were combined for this 
purpose. 

At the central station, No. 2, the amplitudes of the velocity components, 
Us, Ugi.5 and Uj, at the three depths are shown in Table 2. The Table shows 


also the velocity ratios and the time interval by which the current maximum at 
g1°5 and 183 m occurs earlier than at 9m. 


Table 2 
Calibration of North Channel cable: Velocities at central station. 


Us Uo1.5 Uiss 
Date cm/s cm/s cm/s 


Us1.5 


Uiss Time difference 


Us Us Usi.s—Us Uiss —Us 


10 July 165 144 126 0°87 0°77 18 21 
25 July 107 104 95 0°97 0°89 6 33 
Mean 0-92 0°83 12 27 


At station 3, mean depth 130m, it was found that U4¢/U, = 0-86, Ug1.5/U9 = 
0-75, while the phase differences were 7 min and 20 min respectively. 
Considering the variation in depth-mean velocity across the section, 
Ui/U = 0-94, U2/U = 1-05, U3/U = 0°97 
in the above notation. The phase of the current was 15 min earlier at station 1, 
and 3 min later at station 3, compared with_station 2. 

The mean current across the section, U, was compared with the mean cable 
voltage V, as a function of time relative to HW for the same period, July 9, 10 
and 11, as the drift observations. Before this could be done, however, the cali- 
bration of the recording milliammeter in terms of the voltage across the cable was 
required, and this was done simply by introducing a known e.m.f. into the circuit 


containing the cable resistance and the recorder. The calibration factor derived in 
this way was 


1 volt = 149cm/s (2-gokt). 


For July 25, drift observations were available for the central station only, but 
mean values for the whole section were inferred, assuming the same ratio as found 
during the previous period. The calibration factor thus deduced was 


1 volt = 122cm/s (2-37kt). 


ai 

min min 

€ 


The flow of water through the Irish Sea and its relation to wind 273 


Taking the average of these two calibration factors, the mean calibration factor 
is 


1 volt = 135 cm/s (2°62 kt). 


This value may be compared with that found from the Cook measurements. 
From the drogue measurements, the mean current across the section was 0-87 


times the current at 5fm (gm) in the centre of the channel. Applying this factor 
to the Cook observations 


1 volt = 0-87 x 167cm/s = 145 cm/s (2°81 kt) 
which is in fair agreement with the drogue measurements. 

The calibration factor for the St. Margaret’s Bay to Sangatte cable was shown 
to vary from time to time and as might be expected it fluctuated with the seasonal 
variation of the electrical conductivity of the sea water. For the Straits of Dover 
the variation of the calibration factor was investigated in the following way: if 


Ry = range of velocity, Re = range of voltage and Ry = range of elevation which 
was obtained from the tidal predictions for Dover, then 


Having determined Ry/Ry it was then assumed constant, and Ry/Rg was examined 
to determine changes of f. In the case of the North Channel cable the same method 
was adopted, but Ry was found from the tidal levels at Belfast, the nearest standard 
port. Unlike Dover, however, Belfast is not directly on the coast, but is at the head 
of a shallow lough, and as a result the tidal range there does not seem to bear a 
constant relation to the range of currents in the North Channel. This is further 
confirmed by the fact that the harmonic constants for the Belfast tidal elevations 
differ from those of the North Channel currents. For this reason the results 
obtained were not consistent and it was decided to use the calibration factor de- 
rived in 1956 July for all the analysis. The change in calibration expected, due to 


the change in xo/«;, is about 10 per cent, which is no greater than the errors 
involved in the direct calibration. 


Calibration of the Nevin-Howth cable records. No current measurements 
were made during the period of cable signal recording so that the calibration is 
indirect and relies on data tabulated in Tidal Streams, Part II, together with other 
indications of spring rates on Admiralty charts. These estimates are for surface 
flows and were multiplied by 0-87 to give the depth mean velocities. ‘The mean 
current range across the whole section obtained from these values was compared 


with the cable voltage range, averaged over springs for six months, to give a cali- 
bration factor f, of 1 volt = 106-2cm/s. 


5- Effective electrical conductivity of the sea bed 


If the sea bed were perfectly insulating, the potential difference across the 
channel would be equal to the e.m.f. E given by 


E = VH,Lx10%V 


where V = mean velocity across the section, H, = vertical component of the 
Earth’s magnetic field and L = width of channel. The ratio 6 = e/E may then 
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be determined where ¢ is the observed P.D. across the cable. If the sea bed were a 
perfect conductor the cable signal would be effectively short circuited and then 
& = 0, as e = o. On the other hand if no return current could pass through the 
sea bed then 5 = 1, and e = E. Values of 5 for the Irish Sea cables are given in 
Table 3 and show that the shunting effect of the sea bed has reduced the signal E 
by about a half. 

If the cross-section of the channel is represented approximately by a semi-ellipse 
of major axis L and semi-minor axis D, and if xo, «; are the specific conductivities 
of the sea bed and the sea water respectively, then, as shown by Longuet-Higgins 
(1949) 


Values of x; were determined for each area from observations of the temperature 
and salinity of the sea water, and enabled the estimates of xo, given in Table 3 
to be made. 

It will be noticed that although the effective conductivity of the sea bed below 
St. George’s Channel (8-51 x 10-5 [ohm.cm]-') is similar to that below the 
southern North Sea (13-39 x 107-5) and the Straits of Dover (9-74 x 107-5), the 
effective conductivity of the sea bed below the North Channel is approximately 
three times that in the other areas. 


6. Methods of analysis of the records 
(a) General 

The deflection of the recorder traces from the zero line was read off at hourly 
intervals and converted, by a calibration factor determined as in Section 4, to the 
mean water current across the section. By multiplying the current by the cross- 
sectional area, the rate of transport of volume of water through the section was 
found. The long period mean of these values was used as a zero reference in the 
graphical representation of the results. It must be emphasized, however, that this 
is not intended to imply that the drift current was zero when averaged over a long 
period, but is merely a convenient method of representing the variations of the 
flow. 


(b) Tidal currents 

(i) A harmonic analysis was used on a 29 day period of observations on the 
Portpatrick to Donaghadee cable record. The method employed is detailed in the 
Admiralty Tide Tables, Part III (Doodson & Warburg). 

(ii) Periodogram analysis was used on three records of hourly readings over 
29 days. The computations were carried out on a digital computer (DEUCE) 
and for the 696 hourly observations of each series 348 periods between 2 and 696 
hours were tested. The termination of the programme at a trial period of 2 hours 
is a consequence of the observations being taken at hourly intervals, and shorter 
periods could be tested by taking readings at closer intervals. For tidal purposes, 
however, the two hour limit is adequate. With only 696 observations in the time 
series, the resolution of the periodogram is particularly poor in the region of the 
diurnal period, where adjacent test periods are separated by about 50 min. As the 
higher frequencies are tested, however, the resolution improves so that at 2 hours 
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the separation is only 0-3 min. This means that higher harmonics are easily 
distinguishable from the background and there is little interference between the 
constituents. 


(c) Residual currents 


For most of the analysis the main tidal influence was eliminated from the record, 
and an estimate of the daily residual drift was obtained by taking the average of 25 
consecutive hourly values of the current centred at noon. 

Where a detailed analysis for a selected period was required when a high rate of 
flow was observed for a short time, the above method was inadequate and the 
hourly residual current was determined by storm surge analysis (Rossiter 1959). 


(d) Subsidiary data 


One of the main objects of the investigation was to examine the relationship 
between the variations of the flow of water through the North Channel and the 
variations of the sea surface wind. The daily mean winds were required therefore 
and were obtained from the Meteorological Office Daily Weather Reports, where 
wind velocities are given for numerous stations at 6 hourly intervals. For some time 
observations of the wind at West Freugh, a position fairly close to the North 
Channel, were used. Later on, however, measurements at certain times were 
discontinued, and it was decided to use wind data from another source. 

An attempt was then made to compute the sea surface wind from atmospheric 
pressure observations at Malin, Renfrew and Valley. The results obtained, 
however, were higher than any observed winds and it was concluded that the 
method was unsuitable in the present application, partly because of the nearness 
of the stations and partly because of unsteady atmospheric conditions. As an 
alternative the daily means of the winds recorded every six hours at Ronaldsway 
and Prestwick were adopted for the analysis, it being considered that these would 
represent the winds over the northern Irish Sea and the approaches to the North 
Channel. 

The storm surge analysis of the cable record for the period between 1955 
October 31 and 1955 November 6 required more detailed wind data and this was 
kindly supplied by the Meteorological Office of the Air Ministry in the form of 
hourly recorded winds at West Freugh. 


7. Results 
(a) Tidal analyses and periodograms for selected periods 

A harmonic analysis was carried out on a 29 day period of observations of the 
Portpatrick-Donaghadee cable signal between 1955 July 12 and 1955 August 9, 
when the winds over the Irish Sea were fairly light. The method employed is 
detailed in the Admiralty Tide Tables, Part III, and the results in terms of 
amplitudes in cm/s and phases in degrees of the tidal constituents M2, S2, No, 
Ky, O;, Pi, Ma, MS, are given in Table 4. It will be noticed that a large proportion 
of the tide was due to the lunar semidiurnal constituent Ms, which was nearly 
three times larger than its nearest rival. 

Although the above analysis indicates the contributions from the major tidal 
constituents it gives no indication of other oscillations present, whether tidal or 
otherwise, and it was decided, therefore, to carry out a periodogram analysis on 
the same data. The resulting periodogram is shown in Figure 7(a). The pre- 
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dominance of semidiurnal oscillations is immediately obvious as a series of high 
peaks in the region of 12 hours. Diurnal harmonics were also present although 
interference between the individual diurnal constituents was inevitable, due to the 
low resolving power of the periodogram in that part of the spectrum. This poor 
resolution could be improved considerably by taking one year’s observations instead 
of 29 days. A number of other peaks are in evidence and all of them could be 
accounted for in terms of oscillations of tidal origin. Tidal currents are often richer 
in the higher harmonics than tidal elevations, so it was not surprising to find that a 
number of these oscillations, with periods at or near 3, 4, 6 and 8 hours, were 
definitely present in relatively large amplitudes. 

The amount of energy associated with the periods present may be estimated 
from the periodogram by a process described by Davis (1941). When applied to 
the peaks of the semidiurnal tide this method revealed that about 98 per cent of 
the energy of the time series was contained in the harmonic constituents Me, S2 
and No. 

To compare the constituent amplitudes given by the harmonic analysis and 
the periodogram analysis, the amplitude of test periods in the periodogram nearest 
to the appropriate tidal period were noted and are given in Table 4. The amplitude 
ratio shows the closeness of the two sets of estimates. The diurnal constituents 
may be somewhat different from those of the true water movement because the 
cable also detects Earth currents which are known to have oscillations of diurnal 
period. 

The method of periodogram analysis using a digital computer is speedy and 
extensions to the reduction of tidal data in general seem promising.* However, 
the purpose of the present work was to examine the periodicities of the time series 
to decide if there were any periods present which were not of tidal significance. 
In fact all the significant periods appeared to be related to tidal oscillations. 

The same analysis was also applied to hourly observations taken at Portpatrick 
and Nevin during the 29 day period between 1957 July 16 and August 13. The 
results in periodogram form are shown in Figure 7(b) and 7(c), and the amplitudes 
corresponding to common tidal harmonic constituents are given in Table 4. In 
the case of the Portpatrick series the significant periods are present in similar pro- 
portions, so that there were no obviously transient periods in the records. The 
Nevin periodogram shows that the higher tidal harmonics were again present, 
although the diurnal constituents of the tidal currents were more important than 
in the North Channel. 


(b) Daily values of current 

The 25 h mean values of current in cm/s in the North Channel of the Irish Sea, 
referred to the long period mean as a zero reference, are shown in Figure 8(a) and 
(b), where the period of observation extends from 1955 October to 1958 March. 
Similar values of the mean velocity across the Nevin to Howth section (St. George’s 
Channel) are shown (Figure 8(b)), although these measurements were only 
carried out from 1957 January 1 to 1958 March 31. The vertical scale is also 
graduated in units of time rate of transport of water, viz. km3/day, thus taking into 
account the cross-sectional area over which the currents operate. For example, 
with the same mean current of 1ocm/s for the whole cross-section the transport 
across the Nevin—Howth section (effective area 6-879 km?) would be 59-5 km3/day, 


*Mr. D. E. Cartwright of the National Institute of Oceanography has developed the method 
to analyse tidal data in the English Channel. 
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Fic. 7.—Periodograms derived from hourly observations extending over 29 days. 
(a) Portpatrick, 1955 July 12 to 1955 August 9. (b) Portpatrick, 1957 July 16 
to 1957 August 13. (c) Nevin, 1957 July 16 to 1957 August 13. 
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whereas across the Portpatrick-Donaghadee section (effective area 3-933 km?) it 
would be only 33°98 km3/day. 

The magnitude and direction of the long period mean voltage for each 
cable indicated that they were not due to true water movements. For the North 
Channel cable the long period mean would correspond to a current of 2°99 cm/s 
(10-16km/day) in a northwards direction out of the Irish Sea, whereas across the 
Nevin—Howth Cable it would be 7-30cm/s (43-41 km3/day) southwards. These 
results violate the principle of continuity of volume of water because on the 
average, over a long period, water could not leave through both channels. More- 
over, the annual mean northward current through St. George’s Channel, as deter- 
mined from the velocity distribution (Bowden 1950), was 0-35 cm/s, which, as 
well as being an order of magnitude less than the cable estimate, was in the opposite 
direction. It was clear, therefore, that the long period mean voltage on the cable 
could not represent a true water movement, and that no profitable study of the net 
water movements over a considerable period could be made from the cable ob- 


servations. The residual mean voltage on the cables may be attributable to features 
of the Earth current system. 


(c) Correlation of current with regional winds 

The flow through the North Channel is due not only to the local frictional 
effects of the wind on the sea surface but also to the longitudinal surface gradients 
which are produced by wind effects over a much larger area. No computations of 
sea surface gradients have been made in this area as there were no tide gauges 
suitably located. The gradient forces producing flow are mainly secondary wind 
effects and in the following correlation analysis it is the mean wind over a con- 
siderable area of the Irish Sea and its northern approaches which is used, so that 
the effects of surface gradients are taken into account implicitly. 

The method of analysis follows the procedure of Bowden (1956) when a linear 
relation between wind and flow was assumed. The regression equation 


V = cWy+dWe+l (1) 


has two regression coefficients c and d, while V, Wy and Wg are the mean velocity 
in cm/s across the whole section, the northerly wind component and the easterly 
wind component respectively. The last term, J, is that part of the flow which 
appears to be independent of the wind. Equation (1) may be written in the form 


V = AWcos(0—«) +I (2) 


where W is the resultant wind velocity in m/s and A and « are given by the 
conditions: 


= Acosa, d= Asin«. 


From (2) it appears that for a given A, W and «, the maximum value of the mean 
flow V is obtained when @ = a, that is « gives the direction of the wind which is 
most effective in moving water through the North Channel. For the series of 
observations between 1955 October and 1958 March, the constants of the re- 
gression equation are given in Table 5. Under the influence of the most effective 
wind (i.e. from S$ 49° E) the part of the flow dependent on the wind was 0-956 per 
cent of the wind speed. The part of the flow which was independent of the wind 
was apparently 3-36cm/s, but it must be remembered that this included the effects 
of earth currents on the cable, and as a result it is of little value. The compound 
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correlation coefficient Rj,23) gives a measure of the goodness of fit between the 
mean flow and a linear combination of the wind components, and its value in 
this case was fairly high and equal to 0-68. Without doubt there was a significant 
relation between the wind and the water movement. 

Initially it was suspected that there was some inequality in the magnitude of 
the effects of opposite winds. More specifically it looked at first sight as if a 
northerly flow was produced more readily than a southerly one. This might, on 
the one hand, have been due to unequal exposure to all directions at the wind 
measuring stations, or on the other hand, to the effect of the different physical 
characteristics of the northern and southern approach areas of the North Channel. 

To resolve this difficulty the observations were sorted into two groups, in one 
of which the daily mean wind had a component towards N 27°W along the axis 
of the channel, and in the other it had a component in the opposite direction to- 
wards $27°E. The summer and winter periods were considered separately and 
the computed regression equations are shown in Table 5. Although there is a 
difference between seasonal values the two groups are similar, and it was concluded 
that all observations could be given equal weight. 

The question of seasonal variations of correlation also presents itself and to 
investigate this possibility, in the first instance, winter (October-November) and 
summer (April-September) observations were considered separately. The con- 
stants of the regression equations are given in Table 5, and it will be noted that 
the total correlation coefficient Ri23) for the winter series was a good deal higher 
at 0-77 than for the summer series when it was 0-42. Similarly, the value of 
A was somewhat higher at 1-152 in winter than in the summer when it was 0-767. 

As a result the data were further analysed in bi-monthly periods so that the 
correlation variations could be studied in more detail. The computed coefficients 
are given in Table 6, which includes the results of computations from records 
taken between 1954 December and 1955 September, although they have not been 
included in the main analysis. When studying these results it was suspected that 
the recorder biassing arrangements were unsatisfactory and consequently a centre 
zero instrument was fitted in 1955 October to replace the original side zero recorder. 
The immediate effect was to reduce the value of J considerably. ‘The changes of 
the bi-monthly root mean square wind (W2)! and Rj 3) with the seasons are also 
shown. In Figure 9 Ris) is plotted against (W2)! the root mean square wind 
speed, showing that Ris) tends to increase with (W2)!. The seasonal change in 
Ri ,23) is therefore related to the fact that higher wind speeds were experienced in 
the winter than in the summer. 

The values of A shown in Table 6 vary from one bi-monthly period to another 
with extreme values of 0-522 and 1-567 during the main series of observations. It 
was apparent that a high correlation existed between A and (W7?)! and this is illus- 
trated in Figure 1o. A linear regression formula computed under the least squares 
criterion to fit the data was 


A = 0:2313(W2)!—0°386. 4< W2jt< 8. 


The results seem to indicate that it would have been better to correlate the mean 
flow with the square of the wind velocity, to secure a less variable value of A. 
Nevertheless, the method adopted still provides a post prediction flow formula 
where only the wind velocity is utilized, although it should be pointed out that the 
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linear formula tends to overestimate the low flows and underestimate the high 
ones (See also Section 7(d)). 

The values of the coefficient A for the North Channel, in both the summer and 
winter periods, are approximately one third of those for the Straits of Dover. 


Total correlation coefficrent 


5 


Fic. 9.—Variation of the total correlation coefficient, Ri2s) with the 
root mean square wind speed, (W2)#. 


This is probably due to the fact that the North Channel section is 3-3 times as deep 
as that in the Straits of Dover and has a less constricting effect on the flow. The 
ratio of the cross-sectional areas of the Nevin—Howth and North Channel sections, 
from the values in Table 3, is 1-88. In the English Channel, on the other hand, the 
ratio of the cross-sectional area of a section from Portland Bill to Cap de la Hague 
to that of the section in the Straits of Dover is 4°56. 


Wind factor A 


(wt m/s 


Fic. 10.—The wind factor, A, as a function of the root mean square 
wind speed, (W?)!. 


(d) Results of detailed analysis for a selected period 

Heavy gales were experienced over the Irish Sea on 1955 November 2, and the 
effect on the cable signal was easily detected as a displacement of the whole trace 
from its usual mean position. This section describes the detailed analysis of the 
record for the period from October 31 to November 6. The general meteorological 
situation, before and during the period of the storm surge, is interesting, and it 
will now be described briefly before proceeding to the analysis and the results. 
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The sequence of events can be followed on the Meteorological Office weather 
charts which are reproduced in Figure 11 for 0600 on November 1, 2 and 3. 
Although an anticyclone persisted for a few days over the eastern Atlantic before 
October 31 it then moved eastwards across England and a depression off Iceland 
drifted slowly south. On the next day a complex low to the west of the British 
Isles moved south, with some deepening, and was followed by another deepening 
depression which moved south-east from off west Iceland. The wind, which had 
been light over the Irish Sea, increased and became southerly. During the early 
part of November 2 the deep depression moved north, and later on the winds over 
the northern part of the Irish Sea were very strong and at times gusts of over 28-3 m/s 
(55 knots) were recorded at West Freugh. With the movement of the low north- 
wards and then westwards, the winds, which were south and south-easterly, 
moderated, and by next day the gale had abated. It was during November 2, 
when the winds were so strong, that the cable signal was mainly affected. 


Table 6 
Correlation of cable records and mean winds in bi-monthly periods 
Period A a? Rix23) I (W?)t m/s 

*1954 Dec-1955 Jan 1°36 7° 14°32 78 
*1955 Feb—-Mar 1°46 39 0°83 21°72 5°7 
Apr-May 0°35 13 0°36 15°54 
se June-July 0-28 119 0°33 13°83 43 
Aug-Sept 0°47 94 0°39 17"10 48 
Oct-Nov 21 0°69 4°85 47 
Dec—1956 Jan 1°351 61 0-78 6°32 66 

1956 Feb—Mar 0963 47 4°90 6°5 
Apr-May 0°704 77 2°81 4°7 
June-July 0688 39 0°64 2°83 5°5 
Aug-Sept °"759 43 2°21 4°5 
Oct-Nov 32 0°83 
Dec—1957 Jan 1°567 41 0°84 ool 79 

1957 Feb-Mar 0-766 48 0°59 3°55 58 
Apr-May o-611 44 0°45 4°3 
June-July 0°522 47 0°30 6-08 48 
Aug-Sept 0-752 49 0°58 3°16 6:0 
Oct-Nov 1°074 42 0-78 1°74 5°6 
Dec—1958 Jan 1'074 34 2°15 6°6 

1958 Feb—Mar 1:076 33 5°20 6°3 


* Not included in main analysis—-see text. 


The first step in the analysis was to remove the tidal effects from the cable 
record, but where features of flow during a day are required, the rather crude 
method of taking a 25 hour mean is unsuitable. A method of analysis was designed 
by Rossiter (1959) to extract storm surges from tidal records and is applicable to 
the present work. The diurnal, semi-diurnal and quarter-diurnal tides are eti- 
mated and removed, to leave hourly residuals, from which a large proportion of the 
effect of the tide is eliminated. The results for the period from 1955 October 31 
to 1955 November 6 are tabulated and the mean current pulse which occurred 
during November 2 is shown in Figure 12. For the most part the magnitude of the 
residuals was small except during the central day, when aconsiderable non-tidal 
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flow of about 41cm/s, averaged over the whole cross-section of the channel, was 
present. ‘There was also a secondary peak on November 4, which was again closely 
associated with the wind. 

Having obtained the residuals the method of analysis was similar to that 
described in Section 7(c). Results of the correlation analyses of the transport 
and wind over a long period have revealed, that where a linear relation between 
these two quantities is assumed, the total correlation coefficient varies with the 
wind speed (see Figure g), and it is reasonable therefore to correlate the transport 
indicated on the cable with the square of the wind. In these circumstances the three 
variables concerned are: (1) the cable transport, (2) W.Wy, and (3) W.Wg where 
W, Wy and Wg are respectively the wind speed, the northerly component of the 
wind and the easterly component of the wind. 

The computed regression equation was: 


V = 0:0922W? cos(6 + 20° 30’) + 2°058 


where the hourly mean wind velocity W was in m/s and the mean velocity V across 
the section was in cm/s. As in a previous section, it indicates that part of the cable 
signal was not due to wind-induced flow in the North Channel and may be attri- 
buted, to an indeterminate extent, to earth current pick up. The goodness of fit 
of the above formula is shown in Figure 12. 


Cable signal 


T 


i 


00 1000 
1955 November 2 


Fic. 12.—Current surge through the North Channel on 1955 November 
Re cable signal; -----~-- computed from wind data. 


It is interesting to examine the primary and secondary correlation coefficients 
of the foregoing analysis; these are given below: 


riz = 0827, reg = —0'902, 713 = —0°793. 
123 = 0°427, 7231 = —0°719, 713.2 = —O°IQ4. 


The high negative correlation between W.Wy and W.Wg (123.1 = —0°719) 
shows that for the period of observations, wind components towards the north 
were associated with components towards the west, that is, there were persistent 
winds from the south-east quadrant. The low partial coefficient (713.2 = —0°194) 
showed that correlations of transport and W.W, were poor so that the apparently 
high negative correlation given by the primary coefficient (r13 = —0-793) was 
due mainly to the high correlation of the wind components. 
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For comparison the bi-monthly period 1955 October-November gave the 
following correlation coefficients. : 


rig = 0646, reg = —0°361, = —0-437. 
ri23 = 0°582, = —O°1I5, rise = —0°286. 


Here the daily mean wind components Wy and Wz, and not W.Wy and W.Wz, 
were correlated with the water uransport. The wind was much less persistent 
(ro3.1 = —0°115) but, as before, the highest correlation was associated with 
Wy (riz-3 = 0-582). The compound or total correlation coefficient which measures 
the goodness of fit between the transport and a linear combination of the wind 
components was 0°69, which is less than for the quadratic formula where it was 
0°83. 

Referring again to the regression formulae, one would expect a correspondence 
between the predictions of the formula in terms of the wind velocity and the other 
in terms of the square of the wind velocity. In Table 7 the computed daily mean 
flows and the residual movement given by the cable are shown. Whereas in the 
linear formula the wind velocity is a simple daily mean from six-hourly observa- 
tions, in the quadratic formula the mean square of the hourly wind velocities, 
and not the square of the daily mean wind, must be used. 


Table 7 
Comparison of results from prediction formulae 
Daily mean transport velocity across whole of section 


From storm surge 


2shr mean 


bi-monthly linear quadratic 


6-02 


Nov 1 5°04 5°75 10°91 3°45 
2 21°37 21°78 13°84 20°43 
3 5°97 6-41 10°28 7°24 
4 9°37 10°02 10°81 4°80 
5 9°35 


The tabulated values of the daily residual transport velocities determined by 
the storm surge analysis and by a simple mean of 25 hourly observations from the 
cable records, agree reasonably well. The linear regression formula tends to reduce 
extreme values so that low flows are overestimated, whereas rapid flows are under- 
estimated. In the case of the quadratic formula, however, a better fit is obtained. 
During this selected period there was no significant alteration of the daily 
mean sea level at Liverpool. This was due to the fact that S and SE winds do not 
appreciably change the sea level at Liverpool, which responds most effectively to 
a westerly wind acting over the shallow water of Liverpool Bay (Doodson 1924). 


(e) Time dependence of transport generation 

The possibility of the transport lagging the wind effects, although not im- 
mediately obvious from the storm surge data, was also investigated. Total correla- 
tion coefficients were computed for different lag and lead periods between the 
square of the wind velocity and the water transport indicated by the cable signal 
for the same period, 1955 October 31 to November 6. The computations were 
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carried out on a digital computer (DEUCE) for lead and lag periods of 12 hourly 
differences on each side of coincidence. These results are shown graphically in 
Figure 13, and illustrate one rather interesting feature. The highest correlation 
occurs when the transport is correlated with the wind blowing between one and 
two hours previously, and the correlation under these circumstances is 0-85. 
The correlation coefficient differs from that of 0-83 obtained for zero lag by only a 
small amount, but nevertheless the individual cross correlations seem to indicate 


i 
Trausper? leading 


Fic. 13.—Time dependence of the wind produced current between 1951 
October 31 and November 6. 


that the correlation is significantly increased by allowing the wind to lead the 
transport. On either side of the maximum value the correlation diminishes with 
increased time as would be expected. The use of a regression equation associated 
with a two hour transport lag, does not, of course, improve to any marked degree 
the fit of the transport computed by this equation with the observed cable signal. 
It is, however, striking that the water transport through the North Channel reacts 


so quickly to a wind over the North Irish Sea, the lag being not more than two 
hours. 


(f) Correlation of Portpatrick-Donaghadee and Nevin—Howth cable records 


The correlation between the recorded daily mean water movements through the 
North Channel and the St George’s Channel were examined for the period from 
1957 January to 1958 March. For the whole of that time the correlation coefficient 
was 0°55 and by Student's ¢ test it is highly significant. 

The daily values are shown in Figure 8, where corresponding peaks in the flow 
across the two sections are easily discerned. There are, however, long period 
trends associated with the Nevin—Howth record which have no counterpart on 
the North Channel record. These trends are presumed to represent changes in 
the cable zero, as they are too large to correspond to actual transports of water. 
A current of 10cm/s, for example, flowing for 24 hours across the Nevin-Howth 
section, with no flow through the North Channel, would produce a change of 
27m in the level of the whole of the northern part of the Irish Sea. Such 
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trends do not occu= to the same extent in the North Channel record, although an 
apparent northward flow in 1957 June and July is probably spurious, as it had 
no counterpart on the Nevin—Howth record. 

To reduce the effects of long period trends, correlation coefficients for 
3-monthly periods were computed and are given in Table 8. It will be noticed 
that in the summer months there was no significant correlation between the flow 
across the two sections. The winds are, in general, lighter and less persistent at 
this time of year, and the variations in current, particularly across the Nevin-Howth 
cable, were evidently not large enough to stand out above the background fluctua- 
tions. The comparatively high correlations during the winter months, on the 
other hand, confirm the occurrence of periods of wind-induced flow through the 
Irish Sea as a whole. The cable records cannot provide any evidence as to whether 
the flow from St. George’s Channel to the North Channel (or in the reverse 
direction) is concentrated on the western or eastern side of the Isle of Man. 


Table 8 
Correlation of Portpatrick-Donaghadee and Nevin—Howth cable records 


Date No. of observations r significant 

1957 Jan—Mar 66 o'72 yes 
Apr—June 60 no 
July-Sept 74 0°27 no 

Oct-Dec 75 0°49 yes 

1958 Jan—Mar 58 0°78 yes 

1957 Jan-1958 June 333 yes 

(g) Accumulated transport due to wind 

Bi-monthly regression formulae of the form of equation 2, Section 7(b), were - 


given in Table 6 to express the correlation between the wind and the transport 
of water through the North Channel. One term on the right-hand side of this equa- 
tion, namely AW cos(@—«), depends on the daily mean wind vector, while the 
other, J, is not associated with the wind. 


An estimate of the daily transport V’ due to the wind alone may be made 
using the relation 


V’'=V-I 
where V is the daily transport measured on the cable and J is obtained from the 
regression formula for the particular bi-monthly period. 


n=60 
> Vn’ 


n=1 


may then be computed to give the accumulated transport. This has been carried 
out for six bi-monthly periods between 1956 April 1 and 1957 March 31 and the 
results are plotted in Figure 14 where the ordinates of the graph are given in (a) 
km and (b) km per unit area of cross-section. It is revealed that during the year 
a total of 434 km of water flowed out of the Irish Sea through the North Channel 
as a result of the action of the wind, which had an annual mean value of 2:2 m/s 
from the south-west. In terms of a mean flow through the whole section it amounts 
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to 0°35 cm/s or 0-30z2km/day. This current is equivalent to 0-20cm/s through 
St. George’s Channel over the Nevin to Howth cable, and is of the same order of 
magnitude, and in the same direction, as the annual mean flow of 0-35 cm/s 
obtained from a study of the salinity distribution (Bowden 1950). 

This suggests that the mean flow deduced from the salinity distribution may 
be due almost entirely to the residual effect of the wind-induced flow over a long 
period. 

It is also interesting to consider the history of an element of water starting off 
at any instant in the North Channel. It would drift to the south or north depending 
on the wind influence, but the time it took to return to the section would vary 
considerably, depending on the time it started. Figure 13 shows the accumulated 
movement over one year starting from 1956 April 1, and it will be observed that 
the maximum displacements to the south and north were 133 km (1956 December 
14) and 120km (1957 March 23) respectively. The interval between adjacent cross- 
ings of the section may be determined from the same curve by measuring hori- 
zontally along the ordinate corresponding to the day the element left the section, 
to the time when the curve next reached the same accumulated transport value. 
This interval would be very short if the element crossed the section just before the 
transport curve reached a turning point. The interval would be long, however, 
if the element crossed the section at a time such as 1956 April 2, in which case it 
would not have returned from its southward excursion for about 273 days. 

A long displacement in one direction would be expected to result in the element 
of water losing its characteristics, by turbulent mixing with other masses, before 
returning to the North Channel. In this way an irregularly alternating flow could 
result in a “flushing” of Irish Sea water out of the area, even if no net flow were 
apparent. 


8. Conclusions 


The cable method, giving the total flow across a section of the North Channel is 
complementary to other methods of current determination. It provides a continuous 
record of the day to day variations in the flow and enables the correlation of the flow 
with the wind to be established for this area. On the other hand, the cable method 
does not give any evidence on the variation of flow across the section or on the long- 
period mean flow. The variation in current across the section can be calculated 
by Helland—Hansen’s method from the distribution of temperature and salinity, 
but this alone does not enable absolute velocities to be determined. Relative 
currents have, in fact, been calculated from a number of hydrographic sections 
across the North Channel, but the results are not discussed in this paper. Neither 
method is capable of giving a quantitative estimate of the average drift through the 
North Channel over a long period, and for this one must still refer to the estimate 
based on the general distribution of salinity and the influx of fresh water into the 
Irish Sea. 
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Summary 

Theoretical seismograms of Rayleigh waves based on several models 
of mantle structure are compared with actual records for various paths. 
It is found that the model 8099 of Dorman, Ewing and Oliver explains 
seismograms for Pacific paths but does not agree with records from 
Indian—Atlantic ocean paths in the period range shorter than about 
100s. The velocity of the Airy phase corresponding to the group 
velocity maximum is about o-10km/s lower for the Indian—Atlantic 
path than for the Pacific. This difference can be accounted for by re- 
ducing the shear velocity at the top of the mantle under the Indian and 
Atlantic oceans by about 0-1-0-2km/s. The difference between the 
Pacific mantle and the Continental mantle can be explained either by a 
reduction in shear velocity of the low-velocity layer under the Pacific 
ocean or by making the low-velocity zone shallower. 


1. Introduction 


Recently, a comparison of the upper mantle structure under oceans and con- 
tinents was made by Dorman, Ewing & Oliver (1960) based upon the group 
velocity of mantle Rayleigh waves. The existence of Gutenberg’s low-velocity 
layer under continents and oceans was confirmed. However, a significant difference 
in the mantle structure of the Pacific ocean and continents was found. These in- 
vestigators proposed a model for the sub-Pacific mantle in which the top of the 
low-velocity layer is located at a shallower depth than in Lehmann’s model for the 
mantle under continents. Since the shallower low-velocity layer may indicate 
shallower isotherms under the ocean and may consequently explain the unexpec- 
tedly large terrestrial heat flow through the ocean basin, it is important to study 
Rayleigh waves in more detail to compare the upper mantle structure under 
different oceans and continents. 

The upper mantle from the Mohorovitié discontinuity to the bottom of the low- 
velocity layer controls Rayleigh wave dispersion in the range of period from about 
308 to about 200s. Since there is a group velocity maximum in the dispersion curve 
of Rayleigh waves in the above period range (Ewing & Press 1956), the actual record 
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starts with an Airy phase, which makes an accurate dispersion analysis difficult. A 
Fourier method such as applied by Saté (1958) to G waves may yield an accurate 
result. In this paper, we use an alternative method. We start with a layered model 
of the crust-mantle system, compute the dispersion curve for the model, compute 
the impulse response seismogram for the dispersion curve, and compare the com- 
puted seismogram directly with the actual seismogram. We do not seek phase for 
phase agreement as much as signal duration, frequency and arrival time correlation. 
Therefore this is essentially a “‘group velocity method’”’ with the added advantage 
of including the Airy phase in the analysis. Eventually when the initial phase of 
the source and the effect of curvature can be simultaneously allowed for, phase 
and group agreement between the experimental and theoretical seismograms will 
reveal the structure with a greater degree of precision and uniqueness. The disper- 
sion curves are computed by Haskell’s matrix method programmed for the IBM-704 
computer. The computing time is about 3s per layer per period point. The 
impulse response seismograms are computed by the Bendix G-15D electronic 
computer at the Seismological Laboratory, Pasadena. 


2. Rayleigh waves across the Indian-Atlantic oceans 


Dorman, & others (1960) proposed a model named 8099 for the mantle under 
the Pacific ocean based on the group velocity of Rayleigh waves obtained by Sutton 
& others (1960). Aki (1960b) used this model in a study of the mechanism of 
many circum-Pacific earthquakes, and showed that the model explains very well 
the actual features of Rayleigh waves propagated through the North and South 
Pacific oceans. Figure 1 is reproduced from his paper, and shows a comparison of 


No 43 4=42° Aleutian No 18 4=63° South-east Pacific Ocean 


h 


No 52 A=77-5° Japan 4=94° Bismarck Sea 


A= actual T= theoretical - 


Fic. 1.—Theoretical seismograms of Rayleigh waves based on model 
8099 compared with corresponding actual records for Pacific ocean paths. 


actual seismograms recorded at Pasadena and the corresponding theoretical 
seismograms based upon the model 8099. Since he chose earthquakes for which 
the phase angle of the source function is 0 or 7, the observed and theoretical 
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seismograms agree very well, except that the short period waves riding on the long 
period waves of the actual records are cut off in the theoretical seismogram. 

We made a similar comparison for Rayleigh waves propagated through the 
Indian and Atlantic oceans, in order to see if model 8099 also applies to these areas. 
Since we mostly used Pasadena records, the wave paths include continental por- 
tions. For these portions, we used a phase velocity curve that is a combination of 
Press’s curve (1960) and the curve of Dorman & others for the Lehmann mantle; 
the curve is very close to that of model 6£GH which will be described later. 


The impulse response seismograms are computed according to the following 
formula (Aki 1960a), 


fit) = — de, (1) 


where 


A is the epicentral distance, 

w is the angular frequency, 

¢(w) is the phase velocity without the curvature correction, 
din is the instrumental phase delay, 


wy = = is the higher cut-off frequency with 7; = 35s, and 
1 


w = 27 is the lower cut-off frequency with Tz = 200s. 
2 


A comparison of the actual and theoretical seismograms for the Indian and 
Atlantic ocean paths is shown in Figure 2 where the wave forms are arranged not 
according to the absolute arrival time, but in such a way that the agreement is best 
for the longer period waves. We do this because the absolute arrival time of the 
theoretical waves are uncertain owing to a possible error in epicentre. The wave 
group shape, however, is much less susceptible to error in the epicentral distance. 
We see that the observed Airy phase arrives considerably later (by an amount Af) 
than the theoretical one. We cannot exclude the possibility that the observed Airy 
phase is propagated at the theoretically expected velocity, while the observed long 
period waves arrived earlier than theoretically expected. Since this possibility 
seems rather unlikely, we shall assume hereafter that it was not the case. 

We studied Rayleigh waves of 7 earthquakes listed in Table 1, for which the 
great circle paths from the epicentre to the station are shown in Figure 3. The 
delay in the arrival time of the observed Airy phase relative to the theoretical one 
is plotted against the travel distance within the Indian—Atlantic oceans in Figure 4. 
As shown there, the delay time is approximately proportional to the distance 
within the Indian—Atlantic oceans. This fact indicates that the delay of the Airy 
phase is actually caused during propagation through the Indian—Atlantic oceans. 
The delay time is about 6-1 s per ro0ookm. This corresponds to a velocity of the 
observed Airy phase which is lower than the theoretical one based upon model 
8099 by about o-10km/s. 

The next step is to modify model 8099 to account for the observed low velocity 
of the Airy phase propagated through the Indian—Atlantic oceans. 
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New Guinea to Pasadena 


4f 


Indian to Pasadena 


of 
T: Theoretical record (Continent 8099) 
A: Actual record (smoothed) 


Fic. 2.—Theoretical seismograms of Rayleigh waves based on model 
8099 compared with corresponding actual records for Indian—Atlantic 
ocean paths; At is discrepancy in Airy phase time. 


Fic. 3.—Great circle paths across the Indian—Atlantic ocean. 
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Table x 


List of earthquakes 


Epicentre Origin time Region 
(G.C.T.) 
h m s 


40S, 45°5E 13 Indian ocean 
N, 30W 41 Mid-Atlantic 
o'5S, New Guinea 
33S, 179 W South of Kermadec 
36S, 78E Indian ocean 
36°5S, 78-5E Indian ocean 
45, 127°5E Ceram Island 


. Mid-Afiantic to Pasa., Ri 


- Indian Ocean to Pasa.,Ri 


. Kermadec to Lwiro, Ri 
Indian Ocean to Pasa.,R2 


. Ceram to Pasa. R2 
. New Guinea to Pasa..R2 
90 180 
Ain degrees within Indian-Atlantic Ocean 


1 
2 
3 
+ 
5 
6 
7 


Fic. 4.—Delay in arrival time of observed Airy phase with respect to 
theoretical time based on the model 8099. 


3- Models of the upper mantle under oceans and continents 


Before discussing the upper mantle structure under the Indian—Atlantic oc 
we shall enumerate here several mantle models relevant to the discussion. 


1. Lehmann’s model (L). This is based on an unpublished solution by Miss 
Lehmann, derived from her S wave travel-times (Lehmann 1955) for data of 
stations within continents. The dispersion curve of Rayleigh waves for this model 


was computed by Dorman & others (1960). The shear velocity distribution for 
this model is shown in Figure 5. 


2. Gutenberg’s model (G). This model is based on Gutenberg’s study of S 
wave travel-times and is adopted in the Rayleigh wave study by Dorman & others. 
The shear velocity distribution and the dispersion curve of Rayleigh waves for 
this model are shown in Figures 5 and 6 respectively. 
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3. 6EGH. This is based upon the latest result on the shear velocity distribu- 
tion obtained by Gutenberg (1959) from the slope at the inflection points of the S 
wave travel-time curves for different hypocentral depths. This model is identical 
with Press’s model (1960) 6EG for depths smaller than 200km. Case 6EG very 


well explains the observed Rayleigh wave group velocity for Africa for periods 10 
to 70s. 


4. 8099. This is used by Dorman & others to explain the observed group 
velocity of Rayleigh waves for Pacific ocean paths. The model is obtained by 
modifying Lehmann’s model. The shear velocity distribution and the dispersion 
curve for this model are shown in Figures 7 and 8 respectively. 


5. 6EGHP. This model consists of Gutenberg’s mantle (6EGH) and an 
oceanic crust (Figures 7 and 8). 


6. 6EGHP1. This is obtained by modifying case 6EGHP in such a way that 
the shear velocity at the low-velocity layer takes the same value as that in 8099. 
The phase and group velocity curves for this model are almost identical with those 
for 8099 as shown in Figure ro. 


The last three models are related to the mantle under oceans. Cases 8099 and 
6EGHP1 are almost identical and account for the observed Rayleigh wave group 
velocity for the Pacific Ocean. The model 6EGHP, which consists of an oceanic 
crust and Gutenberg’s mantle, gives a Rayleigh wave dispersion curve which 
deviates significantly from the curve for 8099, indicating that Gutenberg’s mantle 
model does not apply to the Pacific ocean basins. 

The dispersion tables computed for this paper appear in Figures 12, 13 and 14. 


4- Mantle under the Indian—Atlantic oceans 


The observed low-velocity of the Airy phase for Indian—Atlantic ocean paths 
can be explained by a modification of 8099 in which either the shear velocity at 
the Moho is reduced or the top of the low-velocity layer is raised. In the model 
6EGHPr’, the shear velocity at the Moho is reduced by 0-1 km/s as shown in 
Figure 9. This reduction gives an Airy phase velocity reduced by about 
0-05 km/s as shown in Figure 10. For case 8099 mod. 3, the top of the low 
velocity layer was raised by 19km compared to 8099, as shown in Figure 9. This 
modification reduced the Airy phase velocity by about 0-07km/s, but also 
lowered the group velocity for periods of 50-125s as shown in Figure 10. From 
these cases, and from many others which we have computed, we conclude that the 
observations can be explained by: (1) assuming lower shear velocity at depths 
from the Moho to about 50km by an amount between o-1 to o-2km/s under 
these oceans than under the Pacific; (2) a combination of (1) together with raising 
the upper limit of the low-velocity zone. These explanations are consistent with 
those of Kovach & Press (1960) who studied Rayleigh waves from epicentres in 
the Indian ocean recorded at the Wilkes station in Antarctica. 


5. Use of the wave shape in deducing fine differences in upper mantle 
structure 


So far, we have been concerned primarily with the group arrival time (e.g. the 


arrival time of the Airy phase) in discussing differences in upper mantle structure. 
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Fic. 5.—Various models of the shear velocity distribution under con- 


tinents. 
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Period in seconds 


Fic. 6.—Theoretical phase and group velocity of Rayleigh waves for 
models shown in Fig. 5. 
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Fic. 7.—Models of shear velocity distribution under oceans. 
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Fic. 8.—Theoretical phase and group velocity of Rayleigh waves for the 
models shown in Fig. 7. 


3 
Upper mantle structure under oceans and continents from Rayleigh waves 299 
‘ | 6EGHP —1_ 
-6EGH 
6EGHP 
| 
|-—-8099 
| 
| 
6EGHP 
\ 8099 | 
CEGHP 
| 8099 | | 
—- - + - + 4+ 
| | | 
| } j 


Keiti Aki and Frank Press 


6EGHP! 


od. 3 — 


\-6EGHP 1 


6EGHPi 


| 
5-0 5.5 


35 =" 
Shear velocity in km/sec 


Fic. 9.—Modified models of shear velocity distribution which account 
for a reduced Airy phase velocity. 
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Fic. 10.—Theoretical group velocity of Rayleigh waves for the models 
shown in Fig. 9. 
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The use of the phase arrival time or wave shape is more effective in deducing fine 
differences in the upper mantle structure. The wave shape, however, depends on 
the space and time factor of the earthquake generating force, and is subject to the 
effect of the curvature of the Earth’s surface as well as to the effect of the “‘pole 
shift” (Brune, Nafe and Alsop, 1961), both of which have negligible effect on the 
group arrival time. Since we cannot precisely specify the earthquake generating 
force at present the use of wave shape for this purpose is premature. It is, how- 
ever, of some interest to see how the theoretical seismograms differ between two 
similar model mantles such as Lehmann’s and Gutenberg’s. 

We shall compute the impulse response f(t) according to Equation 1 for 
Lehmann’s model and for Gutenberg’s model (6EGH) and compare them with 
the Palisades record of a Luzon shock. The shock occurred at 19h 54m 458 G.C.T. 
on 1959 July 18, at Lat. 15°°5 N and Long. 120°-5 E. The great circle path from the 
epicentre to Palisades lies mostly within the continents. A small portion of the 
path is in the Arctic ocean, and we used an oceanic phase velocity for this portion. 
Since the oceanic portion is very small, the choice of any particular mantle model 
for this portion produces a negligible effect on the resultant theoretical seismo- 
grams. We modified case 8099 in such a way that the velocity of the Airy phase 
was reduced by o-10km/s. Figure 11 shows the theoretical seismograms based 


Luzon to Palisades 


—>| imin A= 121-9° 


Lehmann 


Fic. 11.—Comparison of theoretical seismograms for Lehmann’s and 
Gutenberg’s (6EGH) models with Palisades record of Luzon shock. 
Theoretical seismogram is obtained by tracing dots automatically 
plotted by the Bendix computer at the Seismological Laboratory. 


on Lehmann’s model, and Gutenberg’s model (6EGH) and the actual seismogram 
recorded at Palisades. The actual record was smoothed by a symmetric moving 
average, for which the phase shift is zero at any frequency. 

The theoretical seismogram based on Gutenberg’s model shows a good peak 
to peak correspondence with the actual seismogram, while the seismogram based 
upon Lehmann’s model shows a good peak to trough correspondence. ‘This means 
that due to model differences alone, the phase equalization of Rayleigh waves for 
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Fic. 12.—Tabulation of dispersion computations: KD, is dimensionless 
wave number with respect to first layer; D-layer thickness in km, « and 8 
are compressional and shear velocities in km/s; p is density in g/cm‘; 
T is period in s; C and U are phase and group velocity in km/s; 
Mode M1, 1 is Rayleigh mode; N-number of layers used in computation. 
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Fic. 13.—Tabulation of dispersion computations (notation as for Fig. 12). 


6EGHP 6EGHP1 ‘ae 
. 100 
-B40 
- 330 
- 350 
50 “gro 
75 .420 
50 .490 
100 
100 .630 
100 .890 
100 10.170 5.6030 4.130 1¢ . 130 
100 10.585 5.915 - 330 100 330 
100 10.955 6.140 4.490 100 .490 Br 
100 11.275 6. 4.600 10€ .600 
150 11.460 6.38 4.690 150 .690 | 
200 11.755 6.500 4.800 200 
200 12.020 6.610 4.910 200 .910 ig 
S-INF. 12.280 6.740 5.030 3-IN -030 
MODE M1, 1 N= 19 MODE M1, 1 N= 19 
-02605 249.172 4.8400 3.6860 -02596 249.989 4.8400 3.6897 ad 
-O2746 239.346 4.7800 3.6524 4.7800 3.6337 
-02900 229.239 4.7200 3.6350 .02885 230.682 4.7200 3.6128 
219. 4.6600 3.6274 +03050 4.6600 3.6020 
-032 3-600 3-6263 1 +6000 3-3858 
-03471 199. 3.6259 -03437 1.323 3.591 
-O4004 177.512 4.4200 3-239 +0393 180.460 4.4200 3.6070 
-04350 165.639 4.3600 3.6789 256 169.298 4.3600 3.6256 
-04781 152.81 4.3000 3.7101 04641 157 «430 4.3000 3.6487 
-05344 136.638 4.2400 3.7505 .05124 144.604 4.2h00 3.6799 
06138 122.4 4.1800 3.8001 -05762 130.832 4.1800 3.717 
-O7T420 102.772 4.1200 3.8611 .06686 114.0 4.1200 3.7651 
- 10351 74.754 4.0600 3.9447 -08217 94.167 4.0600 
«31524 24.915 4.0000 3.8709 13929 56.385 4.0000 3.97 
10358 74.748 4.0600 "06218 94.163 4.0600 3.8310 
-10834 71.507 4.0550 3.955 .08404 92.185 4.0550 3-838) 
67.810 3.9634 08606 90-135 i «0500 
12187 4.0450 3.9742 - 08826 7.999 4.0450 38935 
-850 4.0400 3.9853 -09068 85.753 4.0400 3.8615 
.14764 pe-737 4.0350 .09 63.392 4.0350 3.8702 
17512 4.516 4.0300 096 4 80.832 4.0300 3.8796 
-21522 36.266 4.0250 3.9991 .09980 78.205 4.0250 3.8900 
24698 31.642 4.0200 3.9734 10373 75.338 4.0200 3-903 
- 26953 29.030 4.0150 3.9472 . 10862 72.037 4.0150 3.914 ae 
- 28684 27.313 4.0100 3- 2263 -11478 68.257 4.0100 3.9282 
- 30214 4.0050 3.8960 12326 6 4.0050 3.9463 

- 31511 24.924 4.0000 3.8712 .13904 4.0000 3.9740 

MODE M1, 1 N= 10 MODE Ml, 1 N= 10 

KD1 T KD1 T U 
-31516 24.921 4.000 3-873 -13904 56.487 4.000 3.9782 
33329 23.612 3.9920 3.82 «31325 25.123 3.9920 3. red 
22.615 3.9840 3.7907 33534 23.515 3.9840 3.8 2 
36197 21.829 3.9760 35342 22.357 3.9760 3.79 
38602 -55 3-9000 3.6671 380 20.826 -9600 3.7032 
. 42289 18.913 3.9280 33063 i8: 5 
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this path would yield source functions which have opposite sense for the two 
models. 

In Figure 11, the seismograms are arranged not according to the absolute time 
but in such a way as to obtain the best fit between the actual and theoretical waves. 
The time shift required to obtain the best fit is indicated by arrows in Figure 11. 
The shift of about 1 minute in Lehmann’s case is probably too great to be accounted 
for by an error in the epicentre location. 


6EGHP1' 
D ALPHA BETA RHO 
5.0 1.520 -000 1.030 
1.0 2.100 1.000 2.100 
5.0 6.410 -700 2.840 
39.0 000 3.3 
25.0 7.840 4.490 3. 
50.0 840) 33.370 
75.0 -020 4.300 3.420 
50.0 8.170 4.440 3.490 
100.0 -500 4.600 3.550 
100.0 9.000 4.950 3.630 
100.0 9.630 2-92 
100.0 10.170 5.630 -130 
100.0 10.585 27923 4.330 
100.0 10.955 -140 490 
100.0 6.288 4.600 
150.0 il. 6.38 4.690 
200.0 11.755 6.500 4.800 
200.0 12,020 6.610 4.910 
S-INF. 12.280 6.740 5.030 
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Fic. 14.—Tabulation of dispersion computations (notation as for Fig. 12). 


This comparison implies at least that Gutenverg’s model 6EGH explains 
Rayleigh wave propagation across the continents as well as Lehmann’s model. 
Since model 8099 for the sub-Pacific mantle was obtained by a modification of 
Lehmann’s mantle, Dorman & others emphasized the difference in the depth of the 
top of the low-velocity layer under the continental and oceanic mantle. If we start 
with Gutenberg’s model (6EGH) in constructing the mantle under the Pacific 
ocean, the required modification (6EGHP1) is not in the depth of the top of the 


| 
MODE M1, 1 N= 19 
-02590 20. 4.8400 4.7155 
02729 4.7800 [f.6310 
-0287 231.2 4.7200 
.03042 222.615 4.6600 
-03232 211.3 4.6000 1.5818 
: -03426 201.952 «5400 
-03657 4.4800 2328 
03923 181.162 4.4200 -6016 
-04237 170.055 4.3600 §§.6197 
-050 145.387 4.2400 §§.6701 
-06597 115.581 4.1200 3.7539 
-08039 .254 4.0600 3.8141 
11721 7.005 4.0000 3.9130 
MODE M1, 1 N= 10 
-11709 67.074 00 3.9136 
4 ‘ 13123 60.0 9920 (3. 9371 
3363 33- 156 9760 3-3372 
a - 27360 28. 238 9680 3.8581 
-29747 26.669 9600 3.8466 
31643 25.122 9520 3.8101 
33269 23.928 9440 3.7706 
-34730 22.982 9360 
.36037 22.194 9280 
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low-velocity zone but in the value of the shear velocity at the lowest velocity layer 
(4°30km/s under the ocean against 4-38km/s under the continent). Thus the 
difference in upper mantle under the Pacific Ocean and continents can be explained 
by a change in velocity in the low-velocity zone as well as by extending the low- 
velocity zone upward. 
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The Correlation of Radio Star Scintillations with 
Geomagnetic Disturbances 


B. H. Briggs 
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Summary 


The correlation between the degree of scintillation of the radio 
source Cassiopeia A and the local magnetic K-index is studied. It is 
shown that the correlation coefficient between these quantities depends 
upon the epoch of the solar cycle at which the observations are made, 
on the position of the source in the sky, and on the radio frequency 
used. For observations on 38 Mc/s made at sunspot minimum, the 
correlation coefficient is small but positive for all positions of the 
source. At sunspot maximum, the correlation coefficient is positive 
when the source is near upper transit and negative when it is near 
lower transit; the average value of the correlation coefficient for all 
positions of the source is approximately zero. It is shown that these 
results can be explained when the finite angular diameter of the radio 
source is taken into account. For some conditions an increase in the 
degree of irregularity of the ionosphere can result in a decrease in the 
degree of scintillation which is observed. It is suggested, therefore, 
that the degree of irregularity of the ionosphere always increases during 
magnetic disturbances in medium latitudes, and that the negative 
values of the correlation coefficient between scintillation index and 
magnetic K-index which are sometimes obtained are due to a secondary 
effect related to the finite size of the source. The results of the present 
paper explain the inconclusive and conflicting results of earlier workers. 


1. Introduction 


Regular observations of the degree of scintillation of the radio source Cassiopeia 
A were made during 1957 and 1958 at Cambridge (52° N, o°W) and Jodrell Bank 
(53°N, 2°W) as part of the programme for the International Geophysical Year. 
The Cambridge observations were made on a frequency of 38 Mc/s and 
those at Jodrell Bank on a frequency of 79 Mc/s. The observations at both 
places were a continuation of earlier work and the methods of observation have 
been described in earlier papers (e.g. Ryle & Hewish 1950; Dagg 1957). The 
records are of the intensity of the source (i.e. the square of the amplitude) and the 
degree of scintillation is conveniently given by expressing the depth of scintillation 
on the record as a fraction of the mean source intensity. At Cambridge an index 
on a scale o to 5 is used to express this quantity, while at Jodrell Bank the results 
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are given as percentages of the mean source intensity; the measurements are of 
essentially the same quantity, and for brevity we shall refer to both as “‘scintilla- 
tion indices’. The scintillations have been shown to be produced by irregularities 
of ionization in the ionosphere, the main effect being due to the F Region. 

It appeared that the observations made on the two frequencies during the 
International Geophysical Year, a period of exceptionally high solar activity, 
would be very suitable for a detailed investigation of the relation of radio star 
scintillations to magnetic disturbances. Previous work on this relationship has been 
inconclusive and to some extent conflicting. Thus Hewish (1952) found no corre- 
lation between the scintillation index observed at Cambridge during 1949-50 and 
the occurrence of magnetic storms. On the other hand Dagg (1957), observing at 
Jodrell Bank during 1955, showed that for a number of individual magnetic dis- 
turbances the scintillation index increased and was quite closely related to changes 
in the horizontal component of the magnetic field. When the magnetic K-index 
was used as an index of disturbance, however, and all the results were taken to- 
gether, the overall correlation was found to be zero. Hartz (1955), observing scin- 
tillations of the same source at Ottawa (45°N, 75°W) and using a frequency of 
54 Mc/s, found only a very small positive correlation between the occurrence of 
scintillations and the magnetic K-index. 

The fact that no clear positive correlation is found between the scintillation 
index and the magnetic K-index is surprising for the following reason. It has been 
shown that there is a close correlation between the occurrence of scintillations and 
the presence of diffuse reflections from the F Region of the ionosphere, when 
observed by the echo-sounding method. This condition is usually referred to as 
“Spread F”’ (e.g. Ryle & Hewish 1950; Little & Maxwell 1951; Briggs 1958a). It 
is therefore reasonable to suppose that the same ionospheric irregularities are re- 
sponsible for both phenomena. Also, it has been shown that in medium latitudes 
there is a strong positive correlation between the occurrence of spread F echoes 
and the magnetic K-index (Shimazaki 1959; Lyon, Skinner & Wright 1960). 
Thus, the degree of irregularity of the ionosphere almost certainly increases during 
magnetic disturbances, and it is difficult to account for the lack of correlation be- 
tween radio star scintillations and magnetic K-index. 

The results of the present investigations are presented in Section 2. It is shown 
that the correlation coefficient between the scintillation index and the magnetic 
K-index is different for the two radio frequencies (38 and 79 Mc/s) and also de- 
pends upon the position of the source in the sky. For observations made on 
38 Mc/s, the correlation coefficient is found to be positive when the source is near 
upper transit, and negative when the source is near lower transit. In view of this 
surprising result, the work was extended to include the sunspot minimum period. 
For this purpose, observations made at Cambridge during 1954-5 were used. It 
is found that for sunspot minimum the correlation coefficient is small but positive 
for all positions of the source. A possible explanation of these results is discussed 
in Section 4, where it is suggested that the scintillation index may be influenced 
by the fact that the radio source has a finite angular diameter. The importance of 
this effect varies with the position of the source and the degree of irregularity of 
the ionosphere, and it is shown that these variations can readily account for the 
fact that the correlation coefficient of the scintillation index with the magnetic 
K-index is sometimes positive and sometimes negative. The same effect would be 
expected to modify the mean diurnal variation of the scintillation index, and this 
point is considered in Section 3. 
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We shall not be concerned in the present paper with the rate of scintillation of 
the radio source, but it may be noted in passing that there is general agreement 
that in medium latitudes the rate increases during periods of magnetic distur- 
bance. This is readily explained, since it is known that the drift velocity of the 
irregularities in the ionosphere increases during magnetic disturbances (Briggs & 
Spencer 1954). Thus, the irregular diffraction pattern formed on the ground is 
earried more rapidly past the observing receiver, and the scintillation rate is in- 
ereased. The high correlation which exists between the scintillation rate and the 
magnetic K-index is probably not significant when the origin of the irregularities 
is in question, though it is of interest in connection with their movements. In con- 
trast, any correlation which may be found between the intensity of the ionospheric 
irregularities and the magnetic K-index is important in the discussion of theories 
of the origin of the irregularities. 


2. The correlation between the scintillation index and the magnetic K- 
index 


In evaluating the correlation coefficient between two quantities like the scin- 
tillation index and the magnetic K-index, it is important to consider the effect of 
periodic variations and trends. Both quantities have diurnal and seasonal varia- 
tiens, and also vary with the sunspot cycle. The diurnal variation alone, which 
is similar for the two quantities with a maximum near midnight, would produce 
@ positive correlation coefficient between the two quantities. This would not 
necessarily imply any close relation between the quantities, since many other 
quantities have a similar diurnal variation (e.g. ground temperature). It is there- 
fore desirable to remove the mean diurnal variations before the correlation co- 
efficient is calculated, or alternatively to calculate the correlation in such a 
way that the diurnal variation has no effect. Similar remarks apply to seasonal, 
sunspot cycle, and any other long period variations which may be present. 

In the present work, these variations were eliminated in the following way. 
Firstly, the data for each month were considered separately. ‘This in itself was 
considered to be sufficient to eliminate seasonal and sunspot cycle variations. 
Secondly, the correlation coefficients were evaluated at fixed times of day by using 
the values for each day of the month. This meant that each value of correlation 
coefficient was calculated from about thirty pairs of values of scintillation index 
and K-index. The correlation determined in this way is between the day to day 
fluctuations of both quantities from their monthly mean values at fixed times of 
day. A disadvantage of this method is that the error of each value of correlation 
coefficient is rather large (about + 0-18), but it has the advantage that the variation 
of the correlation coefficients with solar time and with the hour angle of the source 
can be investigated. 

Magnetic K-indices from the magnetic observatory at Hartland (51°N, 4°W) 
were used. As these are given at three-hourly intervals, it was necessary to work 
with this interval throughout the analysis. The first step, therefore, was to work 
out three-hourly averages of the scintillation indices, which are given at hourly 
intervals for Cambridge and half-hourly intervals for Jodrell Bank. The intervals 
were centred on 01.30, 04.30, 07.30, 10.30, 13.30, 16.30, 19.30 and 22.30 G.M.T. 
The correlation coefficients between the three-hourly averages of the scintillation 
indices and the corresponding K-indices were then calculated from the standard 
formula. In order to investigate the dependence of the correlation coefficient on 
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the hour angle of the source, a mean value for the difference between solar time 
and hour angle was used for each month. This difference is in error by + 1 hr 
at the beginning and end of the month, but this error is not considered to be 
important for the present purpose. Values of the scintillation index were not 
available for every hour owing to equipment limitations; in particular data for 
most of the daylight hours are missing in the Cambridge results owing to the 
serious interference from distant radio stations which is experienced on the fre- 
quency of 38 Mc/s in the day time. However, all values of the hour angle of the 
source are obtained in the course of a year. The data for Jodrell Bank were more 
continuous, with only occasional missing values, but in order to permit a valid 
comparison between the results obtained from the two sets of data, values of the 
index were only used when they were available for both stations. 

The results for 1958 (the later part of the I.G.Y.) are shown in Figure 1(a) 
(Cambridge, 38 Mc/s) and Figure 1(b) (Jodrell Bank, 79 Mc/s). The individual 
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(b) Jodrell Bank , 79 Mc/s. 


Fic. 1.—The values of the correlation coefficient between the scintillation 
index and the magnetic K-index, plotted as functions of solar time, and 
of the hour angle of the source. (a) Cambridge, 38 Mc/s. (b) Jodrell 
Bank, 79 Mc/s. Results for 1958 (mean sunspot number 185). 


monthly values of correlation coefficient are plotted as functions of solar time and 
of the hour angle of the source. Also the mean values of groups of points are shown, 
and these mean values are joined by straight lines. 

Owing to the absence of points during the daytime, no significant conclusions 
can be drawn from the solar plots, though there does seem to be a tendency for 
the correlation coefficient to have a minimum value at midnight. When the results 
are plotted as functions of the hour angle of the source, however, some interesting 
and significant variations become apparent. The clearest variation occurs for the 
38 Mc/s results, for which the correlation coefficient is negative when the source 
is near lower transit, and positive when it is near upper transit; the mean value of 
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the correlation coefficient for all the results taken together is approximately zero. 
The results for 79 Mc/s show a somewhat similar though less marked variation. 
Tike mean value of the correlation coefficient is approximately zero near lower 
transit and has a positive value of about +0-5 near upper transit; the mean value 
far ail the results taken together is + 0-36 + 0-02. 

Im considering the significance of these ‘results, it is necessary to remember 
that the line of sight to the source intersects the ionosphere in a point whose 
lotitude varies with the position of the source in the sky. This point is at its highest 
passibie latitude when the source is at lower transit. For the Cassiopeia source, 
the effective region of the ionosphere for the production of scintillations is quite 
ease to the auroral zone when the source is near lower transit. On these grounds, 
@ higher correlation might have been expected between scintillation index and 
K-index when the source is near lower transit than when it is near upper transit. 
This is exactly the opposite of the variation obtained. Thus, the results cannot be 
explained on the basis of the changing geographical position of the effective region 
of the ionosphere. 

In view of these unexpected results, the analysis was repeated for other years, 
using the Cambridge data for 38 Mc/s only. The results for 1957 and 1959 were 
very similar to those for 1958, and are therefore not reproduced. These were all 
years of high solar activity, near to the peak of the solar cycle. The data available 
for sunspot minimum conditions were incomplete and it was necessary to group 
together a number of months from 1953, 1954, and 1955 in order to obtain sufficient 
data. The results, which refer only to the observations made at Cambridge on 
38 Me/s, are shown in Figure 2. The mean sunspot number for the months which 
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Fic. 2.—The values of the correlation coefficient between the scintillation 

index and the magnetic K-index, plotted as functions of solar time and 

of the hour angle of the source. Cambridge, 38 Mc/s. Results for 1953, 

Aug., Sept., Oct., Dec.; 1954, Jan., Feb., Sept., Dec.; 1955, March, 
June, Aug. (mean sunspot number 13). 


were used was 13. The results may be compared with those shown in Figure 1 
fer 1958, for which year the mean sunspot number was 185. It will be seen that at 
sunspot minimum there is no significant variation of the correlation coefficient 
with the hour angle of the source. All but one of the mean values are positive, 
though the values are very small; the mean value of the correlation coefficient for 
all the results taken together is + 0-10 + 0-02. 

From all the above results, the conclusion is reached that a clear correlation be- 
tween the scintillation index and the magnetic K-index exists only at sunspot 
maximum, and that it varies with the position of the source and with the fre- 
quency used. In particular, for observations on 38 Mc/s, the correlation coefficient 
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is positive when the source is near upper transit and negative when the source is 
near lower transit. It is now clear why earlier workers obtained inconclusive or 
contradictory results since they did not separate the results according to the epoch 
of the solar cycle and the hour angle of the source. Since a positive correlation is 
expected for the reason already discussed, it is reasonable to suppose that the 
negative values obtained when the source is near lower transit are in some way 
anomalous, and produced by the intrusion of some secondary effect. 

In Section 4 a suggested explanation of these negative values will be given, 
based upon the idea that they are due to the effect of the finite size of the radio 
source. However we shall first present a related piece of evidence concerning the 
mean diurnal variation of the scintillation index, which also shows an anomaly 
when the source is near lower transit. 


3- The monthly mean diurnal variation of scintillation index 

Previous workers have discussed the diurnal variation of the scir. ation index, 
and drawn attention to the maximum which usually occurs at a time near midnight. 
In Figure 3 the Cambridge results for 38 Mc/s are presented in the form of monthly 
mean curves for 1954 (near sunspot minimum) and for 1958 (near sunspot maxi- 
mum). The results for 1954 show the expected maximum during the night and 
low values by day. The results for 1958 show a new and unexpected variation. 
For Jan., Feb., March and April, the night time peak is large in magnitude and 
has become double, with a well defined minimum near midnight. The curves for 
Oct., Nov., and Dec., are rather similar. The remaining months show a very com- 
plicated variation with high values by day as well as by night. The feature to which 
we wish to draw particular attention is the splitting of the night-time maximum 
into two peaks. It will be noticed that the minimum between these peaks is most 
clearly defined when the time of lower transit of the source is near midnight. 
(Times of lower transit for each month are shown in Figure 3 by the short vertical 
lines). A careful examination of the curves shows that the minimum does not 
occur exactly at midnight or at the time of lower transit, but at a time roughly mid- 
way between these two times. A similar though less pronounced feature in the 
diurnal curves is also noticeable in the Jodrell Bank results for 79 Mc/s, as re- 
ported by Chivers (1960). In Figure 4 the results for the two frequencies for 1958 
March are compared. It will be noticed that the dip in the curve near midnight is 
narrower and less deep on the higher frequency. 

It is difficult to believe that this secondary minimum during the night, clearly 
associated in some way with the time of lower transit of the source, really represents 
a decrease in the effect of the ionosphere on the radio waves. When the source is 
near lower transit, the thickness of the ionosphere along the line of sight has its 
maximum value, and it would be expected that the irregularities in the ionosphere 
would then have their maximum effect. The results suggest that some secondary 
effect is causing the scintillation index to decrease when the disturbing effect of 
the ionosphere becomes very large. If this secondary effect were absent, the 
diurnal curve for 38 Mc/s might be approximately as shown by the dotted line in 
Figure 4 and this line would represent the true disturbing effect of the ionosphere. 
If this is accepted, it follows that between the points marked (a) and (b) an increase 
in the degree of irregularity of the ionosphere causes a decrease of the scintillation 
index. The negative correlation coefficient between the scintillation index and the 
magnetic K-index obtained when the source is near lower transit could then be 
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ate wneertain owing to the small number of observations made at these 
times. 


1954 and 1958. Cambridge, 38 Mc/s. 
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explained if it is supposed that the degree of irregularity of the ionosphere increases 
with increasing K-index. The results for both frequencies (Figure 1) are quali- 
tatively just as would be expected. 


Scintillation indices (arbitrary scales) 


CMT 
Fic. 4.—Diurnal variations of scintillation index for 1958, March. 
(1) Cambridge, 38 Mc/s. (2) Jodrell Bank, 79 Mc/s. 


4. A suggested explanation of the results: the effect of the finite size of the 
radio source 


The first suggestion that the degree of scintillation of a radio star may be re- 
duced because of its finite size appears to have been made by Little (private com- 
munication). The amplitude of scintillation of the source in Cassiopeia, observed 
at the Geophysical Institute of the University of Alaska, was found to be less than 
for the source in Cygnus, when both were observed at the same angle of elevation. 
A similar result has been reported by Chivers (1960), observing at Jodrell Bank. 
Now the Cassiopeia source is approximately circular in shape with an angular 
diameter of about 4 minutes of arc (Baade & Minkowski 1954). ‘The Cygnus source 
is believed to consist of two small objects each subtending about 1 minute by 0-5 
minutes of arc (Jennison & Latham 1959). Thus, the two sources differ consider- 
ably in solid angle, and if the different parts of the larger Cassiopeia source scin- 
tillate independently, the fluctuations will tend to cancel out. This effect is, of 
course, the radio analogue of the well known optical effect that stars show scin- 
tillations but planets do not. 

In the present work, only the Cassiopeia source is observed. Whether the 
finite size of this source is important or not depends upon the size of the ionospheric 
irregularities and their distance from the observer. We shall now discuss this 
point with a view to explaining the results obtained in Section 2 and Section 3. 

In Figure 5, the ionospheric irregularities of size L are assumed to be con- 
centrated in a layer at a height 4 above the ground. Let ¢; be the root-mean-square 
phase deviation produced across the wavefront of a wave which traverses this layer 
at normal incidence. If the angle of incidence is i, the length of path through the 
irregularities is increased by a factor sec i. Now the root-mean-square phase devia- 
tion is proportional to the square root of the length of the path through the irregular 
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medium (Fejer 1953; Briggs 1958b), and so the phase deviation ¢o for an angle of 
7 is given by 
¢o = di(sec (1) 
Now consider the irregular diffraction pattern produced near the ground on a 
plane normal to the direction of the source (Figure 5). For ¢9 > 1 radian, Hewish 
(1951) has shown that the linear scale of the pattern is of the order of L/¢o. Also, 
fer small angles of diffraction, the pattern turns through the same angle as the 


imetdent radiation. Thus, the diffraction patterns due to opposite edges of the 
tadie source are displaced a distance 20, where 0 is the angular diameter of the 


lonospheric irregularities 


/ Scale Z 


Scale Z/9, 


Fic. 5.—Diagram to illustrate the effect of the finite angular diameter 
of the radio source. 


source and 2 the distance from the irregularities to the observer. For the amplitude 
of scintillation to be significantly reduced due to the finite size of the source this 
displacement must be of the same order as the scale of the pattern, that is 


20 ~ Lido. (2) 
i If the curvature of the earth is neglected, we have sec i = z/h. Thus, using (1) 
: the condition (2) may be written 
dx(sec i)! ~ (3) 


For given values of L, h, and @, the effect of the finite size of the source will be- 
come important when ¢,(sec 7)? is large enough to satisfy equation (3). Thus, the 
effect is likely to be more important when the ionosphere is very irregular (4; 
large), and more important when the source is near lower transit (i = 70°) than 
when it is near upper transit (¢ = 6°). 

The order of magnitude of the quantity L/(4@) will now be estimated. The 
irregularities are known to be elongated (Spencer 1955, Jones 1960), and for the 
present purpose it is their smallest dimension which is relevant. In the work of 
Jones (1960), the sizes of the irregularities were represented by the distance at 
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which an auto-correlation function had fallen to one-half, and a typical value for 
this quantity, measured along the smallest dimension, was 1 km. The values of L 
to be used in equation (3) should be a length which is more characteristic of the 
separation between maxima, and this will be about four times the distance deter- 
mined by Jones. We therefore take L = 4km. For the height 4 a mean value of 
350 km will be assumed. The angular diameter @ of the source Cassiopeia A is 
about 1-2 10-8 radian. With these values, the quantity L/(h@) is found to be 
approximately 10. 

We can now see whether the finite size of the source is likely to be important. 
When the source is near upper transit equation (3) shows that it would be necessary 
for ¢; to be at least 10 radians before the effect would be significant. However, 
when the source is near lower transit, the same equation shows that the necessary 
value of ¢; is then only 2 radians. Previous work (Hewish 1952) has shown that a 
typical value for the root-meai.-square phase deviation on 38 Mc/s for moderately 
disturbed conditions is of the order of 1 radian. For the more disturbed conditions 
encountered during the sunspot maximum year of 1958, values of the order of a 
few radians would be a reasonable assumption. It therefore appears that the 
effect of the finite size of the source is likely to be important on 38 Mc/s only when 
the source is near lower transit. Now the root-mean-square phase deviation is 
inversely proportional to the frequency, and so for 79 Mc/s any effects which may 
be attributed to the finite size of the source are likely to be less pronounced, and 
to be confined even more closely to the time of lower transit. 

In general, the effect will be a maximum when the quantity ¢;(seci)? is a 
maximum. Now 4; reaches a maximum at some time near midnight, but (sec 7)? 
is a maximum at the time of lower transit. Thus the greatest effect will be observed 
when lower transit occurs at midnight, because both terms then reach their 
maximum values simultaneously. If lower transit is not at midnight, the greatest 
effect will be observed at some time between the time of lower transit and mid- 
night, when ¢;(sec 7)? is a maximum. 

To explain the results of Section 2, it is only necessary to note that when the 
source is near lower transit there will be occasions when an increase in the degree 
of irregularity of the ionosphere (i.e. in ¢;) will enable equation (3) to be satisfied 
when previously it was not, and will therefore result in a decrease in the amplitude 
of the scintillation. Thus, if the irregularity of the ionosphere increases with in- 
creasing magnetic K-index, the correlation between scintillation index and K- 
index may, nevertheless, be negative. This effect will be more pronounced on 
38 Mc/s than on 79 Mc/s. The results (Figure 1) show that on 38 Mc/s the corre- 
lation is negative for a considerable period around the time of lower transit. On 
79 Mc/s however, the individual values of the correlation coefficient for the times 
near lower transit are sometimes positive and sometimes negative, with a mean 
value of approximately zero. This is just the type of result which would be 
expected. At sunspot minimum the effect disappears, probably because ¢; is never 
large enough for equation (3) to be satisfied, even near lower transit. 

The results of Section 3 can also be explained by the decrease of scintillation 
index when the finite size of the source becomes important. The main features of 
the curves (Figures 3 and 4) are as would be expected, including the time at which 
the subsidiary minimum occurs. As previously mentioned, this minimum occurs 
between the times of lower transit and midnight, and is most pronounced when 
these times coincide. It is absent at sunspot minimum. 
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5- Further discussion 


It has been shown that the finite size of the Cassiopeia source can explain the 
unexpected results obtained when the correlation of the scintillation index with the 
magnetic K-index is studied. It seems likely that for latitudes near 50°, the degree 
of irregularity of the ionosphere always increases during magnetic disturbances, 
and that the negative values of the correlation coefficient which are sometimes 
obtained when the source is near lower transit are due to this secondary effect 
related to the size of the source. 

It should be mentioned that at low latitudes, within 20° of the equator, radio 
star scintillations and spread F echoes are both negatively correlated with the 
magnetic K-index (Wright, Koster & Skinner 1956; Koster & Wright, 1960; 
Shimazaki 1959a, b; Lyon, Skinner & Wright 1960). Thus, near the equator, the 
degree of irregularity of the ionosphere decreases during magnetic disturbances, 
in contrast to the increased irregularity at higher latitudes. 

It appears from the present results that the scintillation index of the Cassiopeia 
source as measured on 38 Mc/s is not a very useful measure of the degree of irregu- 
larity of the ionosphere. Results obtained on 79 Mc/s are more reliable in this 
respect, but even on this frequency the finite size of the source is sometimes 
important. The degree of irregularity of the ionosphere is probably not very well 
represented by the mean diurnal curves of scintillation index of the type shown in 
Figure 3, and theories of the origin of the irregularities should be not based upon 
their detailed shapes as was, for example, the theory of Harrower (19572, b). 

In an earlier paper, (Briggs 1958a) the correlation between the scintillation 
index and the spread F echoes observed in vertical sounding of the ionosphere was 
studied. It was shown that the highest correlation was obtained when the iono- 
spheric observing station was directly below the irregularities which lay in the 
line of sight to the source. This result was as might have been expected, and it 
was used to make deductions about the sizes of the patches of irregularities. These 
results were for 1954 and 1955, a period near sunspot minimum. When the calcu- 
lations were repeated for a period near sunspot maximum, a completely different 
result was obtained. When the source was near lower transit, the correlation be- 
tween scintillation index and spread-F index was often negative, even when the 
observing station was almost directly beneath the irregularities. It is clear that 
this result can also be explained as due to the finite size of the source. Under such 
conditions, the spread-F index is a more reliable index of the degree of irregularity 
of the ionosphere than is the scintillation index. 
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Preliminary Palaeomagnetic Measurements on Silurian 
Sediments from South Africa 
K. W. Graham and A. L. Hales 


(Received 1961 March 16)* 


Summary 

This paper describes palaeomagnetic measurements made on the 
Lower “Shales” of the Table Mountain Series, near Cape Town. The 
age of this formation is probably Silurian. Samples collected at the 
surface were found to be inconsistently magnetized, but those from 
a fine red mudstene bed exposed in a deep read cutting gave consistent 
results. Bight samples from this bed give a mean direction of 
D = 341°8. I = 3°-5 with the south pole downwards. 


1. Introduction 


Considerable difficulties have been found in interpreting palaeomagnetic 
measurements on surface samples of South African rocks. It is thought that a 
major part of the difficulty in interpretation arises from the remagnetization of 
near surface samples by the magnetic fields associated with lightning currents 
(Graham, submitted for publication). At depth these difficulties do not arise and 
significant results have been obtained. (Gough 1956, Graham & Hales 1957). 

Many of the significant results obtained elsewhere were from samples of red 
sandstones and shale. For this reason an extensive sampling programme was 
carried out on the red sandstones and mudstones on the slopes of Chapman’s 
Peak near Cape ‘Town. 

These sediments are known as The Lower Shales of the Table Mountain Series 
of the Cape System. 


2. The age of the Cape system 
The Palaeozoic “Cape System” of South Africa has the subdivisions: 
Witteberg Series 
Bokkeveld Series 
Table Mountain Series. 
The Witteberg Series contains occasional plant fossils of Devono—Carboniferous 
affinities. The Bokkeveld Series contains abundant marine foss.". of Devonian 
age. Since no identifiable fossils have been found in the Table Mountain Series, its 
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age must be inferred from that of the Bokkeveld Series. Du Toit (1954) regards 


the T.M.S. as “certainly not younger than the Lower Devonian’’, and a Silurian 


3- The sampling 


In the first sampling a number of hand samples and 6-in cores were collected 
from the Lower “Shales” of the Table Mountain Series at Chapman’s Peak (Lat. 
33°°1 S. Long. 18°-3 E.) near Cape Town. The directions of magnetization of 
these samples are shown in Figure 1(a). It was thought that the near-random 
patterns of magnetization of the first set of samples might be due to slight weather- 
ing, or thermo-stress cycling of the outermost few inches of the outcrop and that 
more consistent results might be obtained if deeper samples were obtained. 
Accordingly a portabie drill rig was developed with which oriented cores could 
be obtained to a depth of three feet. (Graham & Keiller 1960). 

Using this rig, 46 3-ft cores were drilled from rock exposed in a deep road 
cutting as well as on the mountain side below the cutting. Some of the cores were 
drilled upwards or downwards at small angles to the horizontal so as to pass very 
slowly through a few inches of the untilted sediments. Other cores were drilled 
at steeper angles to sample two feet or more of the geological column. Figure 2 
shows diagramatically the relative positions of the various boreholes and the 
vertical thickness of sediment represented by each of the 3-ft cores. The portion 
of this diagram above the dotted line represents the road cutting while that below 
the dotted line represents the steep mountain side. The wall of the road cutting 
from which the samples were drilled is about 50 ft to the south east of the natural 
cliff and almost parallel to it. 


4- The samples from the cutting 


Three pairs of cores were drilled about fifty feet apart so as to intersect exactly 
the same few inches of a red siltstone bed about two feet thick. Another core (39) 
was drilled so as to cut the entire bed. Core No. 10 also came from this bed. Other 
cores were drilled from the almost “se ge quartzite exposed in the road cutting 
above and below the fine grained bed 

The directions of magnetization of the samples from the red siltstone in the 
cutting are extremely consistent both laterally over the 10oft sample as well as 
over its thickness. As examples, the directions of magnetization of cores 17, 39 
and 19 are plotted stereographically (on the lower hemisphere) as Figures 1 (b), (c) 
and (d). The mean directions of magnetization of the mudstone cores, are given 
in Table 1. The semi-angles of the cones of 95 per cent confidence (<) 
were calculated by Fisher’s method (Fisher 1953) and are also given in Table 
1. The mean direction of magnetization of this layer, giving each core unit 
weight, is declination 341°-8, inclination 3°-5 with a South pole downwards. 
The semi-angle of the cone of 95 per cent confidence («) is 4°°5. 

The quartzite above and below shows considerable scatter in the directions of 
magnetization of the specimens from each core, but qualitatively the directions are 
not inconsistent with those from the fine grained bed as may be seen from Figure 
1(e) and 1(f). 
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Fic. 1.—Lower hemisphere stereographic plots of the directions of 
magnetization of (a) the reconnaissance samples; (b) to (f) the speci- 
mens from cores 17, 39, 19, 43 and 48 respectively. 


@ = South pole downwards x = North pole downwards. 


gh 
3 
i 
e) 


re! 126 


6 
8 
9 
z 


or 720 
jezuozisoy 


Ayjequozisoy 4,05 Jnoge | 


peoy 
paq pauyesd our 


ie Palaeomagnetic measurements on Silurian sediments 
on from South Africa 321 
7 
OW 
1 
| | 
i 
3 
= 
1 
: 
1 
1 
= 
‘ 
\ 
\ = 
ia 
1 
‘ 


322 UK. W, Graban Bates 
or 
+ 
N 
* 
| 
Kor 
2 
4 Fic. 3.—Lower hemisphere stereagraphic plots the directions sii 
magnetization of the specimens fram cores Nes. 7, 8, 35, 36, 44 and 46 es 
respectively. @ South pele downwards: North pele downwards. 


Palaeomagnetic measurements on Silurian sediments from South Africa 323 


Core No. No. of specimens D (deg) F* (deg) % (deg) 
10 38 341°7 +5°2 4°5 
64 


15 


$'5 
5°7 
5°4 
28 


+02 6-2 
* A positive inclination indicates a South pote downwards. 


5. The samples from the surface outcrop 

Broadly speaking, the cores drilled from the fine grained red, arkosic sandstone 
forming the surface outcrops show two prominent directions of magnetization 
with intermediate directions. The first direction is very similar to that obtained 
from the samples inside the cutting, and consists of almost horizontal South poles 
to the NNE. The second prominent direction consists of North poles nearly 
vertically downwards. 

Samples distributed both horizontally along the same sedimentary layer and 
vertically up the geological column show that these surface outcrops are incon- 
sistently magnetized. In one spectacular instance, samples Nos. 7 and 8, although 
very close to each other both vertically and horizontally (see Figure 2, were 
magnetized in widely differing directions. Cores Nos. 35 and 36 were drilled 
vertically above numbers 7 and 8 respectively, in such a way that each would cover 
the material sampled by both 7 and 8 as well as the vertical interval between them. 

No. 35 agreed in its direction with No. 7 and No. 36 with No. 8. Stereographic 
projections of the directions of magnetization of the specimens from each of cores 
7, 8, 35 and 36 are shown in Figure 3 (a), (b), (c) and (d). Figure 3(e) and (f) 
represents cores 44 and 46 and serves as a further example of horizontal incon- 
sistency in the directions of magnetization of the surface outcrops. 


6. Discussion 


Since the directions of magnetization of the surface outcrop cores from the 
same horizons and only a few feet apart differ so widely, it is clear that some local 
cause has been effective in magnetizing them. It follows that the measured 
directions do not give any meaningful information on the direction of the Earth’s 
field at the time of deposition of the sediments. It is possible that the specimens 
have been remagnetized by lightning currents running down the cliff towards the 
sea. 
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On the other hand, allowing for the éxpetted poor grouping of the directions of 
magnetization of the coarse grained qnartzites, samples collected from a deep road 
cutting are consistently magnetized, secular variation may net have been 
meaned out completely, it seems probable that the mean direetion of magnetization 
of the red, fine grained bed is a reasonable appreximatien to the direction of the 
field during Silurian times: - 

Assuming an axial dipole field, the South pele relative to Africa would have the 
present-day coordinates of longitude 10°-g W, latitude 50°-3N (or the North 
Pole 169°-1 E, 50°-3S). The co-latitude of the Cape Peninsula relative to this 
Silurian pole would be 88°, which places South Africa very nearly on the equator 
at that time. The presence of a tiflite in the Upper T.M.S. Shales is puzzling and 
demands further work on the Cape System. 

It is interesting to compare this pole pesition with these found for other con- 
tinents. Figure 4 is based on Cox and Doell’s Figure 23 (Cox & Doell 1960) and 


Fig. 4.—Early Palaeozoic Virtual Geomagnetic Poles (after Cox and 
- Doell 1960, Figure 23). The numbers and symbols are those used by 


Cox and Doell. 
© European D = Devonian 
( North American S = Sifurian 
Australian = Ordovician 
) Asian @ = South magnetic poles 
© African O = North magnetic poles 


® = Poles from measurements with mixed polarity. 


shows clearly that the pole position found from the Lower Shales of the Table 
Mountain Series is not in agreement with those of Ordovician, Silurian er Devonian 
age from any of the other continents. 
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Infra-red Physics 
| Volume 1, No. 1, 967 
(Pergamon Press. Quarterly, annual subscription £7.) 


Infra-red Physics, published by Pergarbon Press, has given your reviewer ex- 
perience of the razor-edge of ambivalence. There is pleasure that this publicly 
neglected subject has a Journal of its own, herrer that yet another journal must 
be purchased by libraries and:scanned by readers. There is delight in the progress 
of the subject over the last decade; regret at its regressions. 

Wild surmise greeted the ahnauncement by Hanbury-Brown and Twiss of the 
positive result of their famous experiment on photen-correlation. This shows 
how teaching still largely neglécts both physical statistics and infra-red physics. 
The worker in infra-red (which is our present concern) is daily confronted with 
the wave-particle duality, because the infra-red spectrum spans the transition 
from the Rayleigh-Jeans realm of mierowaves to the classical-particle statistics of 
the optical region. It has therefore a special claim to the attention of students of 
fundamental physics. 

Volume 1 of the Journal opens with an histerical paper on Sir William Herschel. 
There is a welcome review article on infra-red detectors. A number of articles 
describe devices; these inelutte sources for communication, tellurium detectors, 
germanium dioxide glasses, interference modulation of millimetre waves, a sun- 
seeker, phophor photography, absorption-edge thermal imaging, and a scanning 
Fabry-Perot. There is a researeh note on absorption in lead halides. 

A notable advance, as compared with the literature of 10 years ago, is the uni- 
versal use of physical units of deteetivity. The days of standardization by cigarette 
ends, soldering irons and signal sehool black have passed unmourned. 

A regression is that particle-statistieal approximations, using crude estimates of 
effective wavelength, have often replaced the more exact formulae for ultimate 
sensitivity. The term “‘ultimate-sensitivity” is due to Jones (1947), and “detectivity” 
was suggested by your reviewer (Jones, tg52). “Ultimate sensitivity and practical 
performance” is a direct quote frem the reviewer (Fellgett 1949), who is disap- 
pointed that the article under this tithe did not use, or refer to, the more exact 
formulae for estimating these quantities, which it is now believed have been cor- 
rectly extended to non-equilibrium conditions (Fellgett, Jones & Twiss, 1959). 
The loss of accuracy is small at present, but the loss of msight is great. 

Science has so far made. relatively little use of the powerful conceptual and 
technical methods now available for the creation, copying, and retrieval of docu- 
ments. Each new journal raises the questien of how long conventional publication 
can continue to be a major method of seientific communication. ‘There are cheaper 
means of creating easily-read text than by printing, and more effective means of 
search than by using subject-indices or seanning the pages of a journal devoted 
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to a particular aspect ef science. Yet the prestige and -‘nity of a published paper 
have an ineeplaceabic valuc. Reforesing contrifutes «ssentially to this value, but 
may make easier te publish a mediocre paper than either a very bad or a very 
goed ene. An author whe fails to place his manuseript with one journal can‘.try 
elsewhere, and in this way can obtain an assessment which is disinterested, and is 
seen to be so. Unified publication is unacceptable because it would destroy this 
freedom. 

As scienc* expands, new journals are needed. But if their number continues‘to 
mulhiply unchecked, the whole system may break down. ‘The abligations of libra- 
ries have the effect of blunting economic selection. Some curb must be found, 


and it must be closely related to the exploitation of new techniques of scientific 
communication. 
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1. Authors are advised to follow the general suggestions “or the preparation of 
scientific papers given in Notes on the preparation of papers to be communicated to 
the Royal Society and in The printing of mathematics by T. W. Chaundy, P. R. 
Barrett and C. Batey (Oxford University Press, 1954). In particular the abbrevia- 
tions for the names of periodicals used in the Geophysical Journal are those given 
in Notes on the preparation of papers . . . and the symbols, signs and abbreviations 
are those recommended for British scientific publications by the Symbols Commit- 
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2. Papers should be typewritten with double specing and on one side of the 
paper only. Two copies, one of them the top copy, should be submitted. All papers 
must be accompanied by a summary; summaries are not required for Letters to 


3. Tables are printed without rules. They should be numbered serially with 
Arabic numerals. Table headings should be brief. Units should be placed at the 
head of the column. Tables should be typed on separate sheets and their positions 
in the text indiceted on the copy. 


4. Illustrations, especially photographs, should be kept to a minimum, and the 
same information should not in general be given in both tables and illustrations, 
Line drawings should be in dense black ink on smooth white board or transparent 
tracing film. Original drawings must not be lettered; lettering should be indicated 
on copies or photoprints which must be supplied in addition to the origimils. 
Wherever possible, lettering should be kept outside the diagram so that it may be 
set in type. Diagrams should be drawn at twice to three times the size at which 
they will be printed. The maximum dimensions for diagrams, when reduced and 
with their lettering and captions, are 8 in by 5 in. A wide variety of flat tinte of 
dots, lines and shadings can be applied to line drawings by the blockmaker. Photo- 
gtaphs for reproduction should be unmounted glossy prints and should be accom - 
panied by lettered prints. 

5. References should be quoted in the form recommended by the Royal Society 
but authors may also quote titles of papers if they wish. References in the text are 
made in the form (Smith 1940). 


6. Authors whose work involves much mathematics should follow carefully the 
recormmendations in The Printing of Mathematics. 
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